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ABSTRACT 

Young circumstellar disks that are still embedded in dense molecular envelopes may differ from their older counterparts, but are 
historically difficult to study because emission from a disk can be confused with envelope or outflow emission. CO fundamental 
emission is a potentially powerful probe of the disk/wind structure within a few AU of young protostars. In this paper, we present 
high spectral (R = 90, 000) and spatial (~ 073) resolution VLr/CRIRES M-band spectra of 18 low-mass young stellar objects (YSOs) 
with dense envelopes in nearby star-froming regions to explore the utility of CO fundamental (Av = 1) 4.6 jim emission as a probe 
of very young disks. CO fundamental emission is detected from 14 of the YSOs in our sample. The emission line profiles show a 
range of strengths and shapes, but can generally be classified into a broad, warm component and a narrow, cool component. The broad 
CO emission is detected more frequently from YSOs with bolometric luminosities of < 15 L Q than those with > 15 L Q . The broad 
emission shares many of the same properties as CO fundamental emission seen from more mature disks around classical T Tauri stars 
(CTTSs) and is similarly attributed to the warm (~ 1000 K) inner AU of the disk. CO emission from the inner disk is not detected from 
most YSOs with a high bolometric luminosity. Instead, the CO emission from those objects is produced in cooler (~ 320 K), narrow 
lines in 12 CO and in rarer isotopologues. From some objects, the narrow lines are blueshifted by up to ~ 10 km s , indicating a slow 
wind or outflow origin. For other sources the lines are located at the systemic velocity of the star and likely arise in the disk. For a 
few YSOs, spatially-extended CO and H2 S(9) emission is detected up to ~ 2" from the central source and is attributed to interactions 
between the wind and surrounding molecular material. Warm CO absorption is detected in the wind of six objects with velocities up 
to 100 km s -1 , often in discrete velocity components. That the wind is partially molecular where it is launched favors ejection in a 
disk wind rather than a coronal or chromospheric wind. 

Key words. Proplanetary disks; Stars: protostars: profiles; infrared: stars; Techniques: spectroscopic 
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1 . INTRODUCTION support in hydrodynamic models of collapsing cold cores dBatel 

OW- ♦ 11 A-v 1 h.1 1 • a a r 2010;Machida&Matsui^|20li. 
j^j Circumstellar disks play a central role in the growt h 01 young 

stars and in the formation of planetary systems (e.g. IShu et all the P rotostar is stin embedded in an optically-thick 

*i ■ [19941 iLissauerlflUh . Disks channel material onto the star and envelope, the disk is difficult to study because emission from 

VO 1 provide an environment for grains to grow into planets. While the dlsk ' lf P resent > can be confused with emission from the 

O ! the star+disk system is still enshrouded in the envelope, the envelope. Observations with high spatial resolution can break 

^ ■ star builds up most of its mass, likely accreted from the disk. thls degeneracy by spatially discriminating between the loca- 

Tl ' Meanwhile, the disk mass can be replenished by the surrounding tlon f compact disk emission and extendgdenyelope emission 

> ! envelope. Eventually the envelope and later the disk disappear, Ig£. Padsett et all 1999] | EisneretaU |2Q05j | Beckford et all 

leaving only the star, any planets or planetesimals that may have EQQS | Tobinetall |2010|). For example, sub-mm interferometry 

formed and remnant debris nas been use d to separate a compact emission source, as ex- 

During the Stage phase of young stellar object (YSO) evo- P ected for a disk ' and an extended emission source, as expected 

lution, the dense envelope surrounds the star and any disk that foran envelope (eg. |Keene & Massonj |199d; |Wilner et al. , 

may be present. Most of the luminosity from the central pro- l996 > Loone y et al " 200 °; Brown et al " 200 °; Jensen et al., 

tostar is reprocessed by col d dust and esca pes as far-IR emis- ™ l2Q07t | Lommenetall |2008fc | j0rgensen et all |2009j). A 

sion (e.g. Andre et al., 1993; Robitaille et al., 2006; Crapsi et al., wlde ran 8 e of disk-to-envelope mass ratios have been mea- 

2008). At the Stage 1 phase of YSO evolution, the envelope sured for ^-mass young stellar objects, as would be expected 

and disk masses are similar. Any disk likely forms during the for a sample that includessystems at a range ofevolutionary 

Stage phase. Eva nT^td] (12009) measured population fractions P hases flLommen et all |2008j | j0rgensen et all |2009j). In most 

and converted them to timescales to estimate that the Stage cases the disk+envelope mass is much smaller than the stel- 

stage typically lasts 0.1-0.16 Myr and that Stage I stage typi- lar mass ' indicating that the main accretion phase has ended, 

cally lasts ~ 0.5 Myr. Fast disk formation also has theoretical Dlsk masses are similar for objects in Stage and Stage I 

stages ( J0rge nsen et all 120091) . which suggests that once the 

* This work is based on observations collected at the European disk grows to a certain size, the mass that the disk loses to 

Southern Observatory Very Large Telescope under program ID 179.C- the star is replaced by accreting mass from the envelope. The 

0151. presence of disks around Stage I objects is supported by obser- 
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vations of Ke plerian emission profile s in cold rotational lines 
of HCO + (e.g. iHogerheiide & Sandelll l2000t LSaito et ail 1200 It 
iBrinch et all l2007t iLommen et aUl2008l) . 

These (sub)-mm observations trace molecular emission from 
the disk at large radii, but little is known about the inner AU 
of young disks that are still embedded in envelopes. If the 
region is similar to the inner disks around classical T Tauri 
stars (CTTSs), the chemical and temperature structure can be 
probed by the same line diagnostics used to study CTTS in- 
ner d isks (see reviews by iNaiita etafl . l2007t iDullemond et all 
120071) . Unfortunately, many of these techniques are difficult to 
apply to embedded objects because emission from the disk, en- 
velope, outflow, and cloud can be confused and because embed- 
ded objects are highly reddened. For example, in Spitzer IRS 
spectral surveys of young stars, emission in mid-IR gas lines 
is detected more frequently and is more lu minous from embed- 
ded objects than from classical T Tauri star s dLahuis et all 2007 ; 
Flaccomio et al l l2009t Lahuis et al .,2010). However, follow-up 
of potential disk tracers, in cluding \Q II, r Ne III and H? emission 
at high spectral resolutio n (Hartigan et al., 1995 ;lHerczeg et all 
2006; Ivan Boekel et all l2009t iPascucci & Sterzikt l2009h. at 
high spatia l resolu tion (I Walter et all 12003: Saucedo et al., 2003; 
Beck et al], 120081) . or in large samples (Giidel et al.. 1201 Ol) . in- 
dicate that even for CTTSs, outflow emission often dominates 
over any disk emission. Emission from disks and outflows can 
be even more confused for Stage 1 objects because they typi- 
cally drive more powerful outflows than CTTSs. Indeed, emis- 
sion in H2 1-0 S(l) and 2-1 S(l) lines from embedded objects 
shows significant spatial extent but with small cen tral velocities , 
indic ating that the gas is kinematically quiescent dGreene et all 
120101) . H2O emission in the near-IR and sub-m m has also been 
attributed to disks around embed ded objects dCarr et all [2004; 
JOrge nsen & van Dishoecktl20loh . 

Fundamental-band (Av = 1) CO emission in the M-band 
has been demonstrated to be a powerful probe of warm (~ 500- 
1000 K) gas in terrestrial planet forming regio ns of disks around 
CTTSs dNaiita et al.Ll2003llRettig et alfeooilGoto et all f2006: 



ISalvk et alll2007l l2009). For a Keplerian disk, the spectral width 
of the emission, when combined with a known disk inclination, 
constrains the lo cation of the inner radius where the CO gas is to 
~ 0.05 - 0. 1 AU dNaiita etalll2003h . near where CTTS disks are 
expected to be truncated by the stellar magnetosphere (fK oni gl , 



\l99^ IShu et all fTf^l Ijohns-KrullL l2007t iGregorv et all 1200 
Spectro-astrometry of CO fundamental emission from three 
CTTS disks with inner dust holes confirms that the CO emis- 
sion is produced in a disk and demonstrates that molecular 
gas is located within the dust hole (Pontoppidan et al. 2008; 
see also Pontoppidan et al. 2011). CO fundamental emis- 
sion i s also commonly detected in disks around H erbig AeBe 
stars dBlake & Boogerttl2004tlBrittain et alll2007l) . The Herbig 
AeBe disks with inner dust holes also lack warm CO emis- 
sion within those dust holes, in contrast to CTTSs (|Goto et all 
2006; Ivan der Plas etafl l2009t iBrittain et all 120091) Several 
Herbig AeBe stars also show bright emission from highly- 
excited vibrational levels, up to v = 7, likely indicative of strong 
UV emission from the s tellar photosphere (Brit tain et all 12009: 
Ivan der Plas et all 120091) 

Several studies have shown that warm CO is also present 
around embedded objects. In one case, GSS 30, CO fundamen- 
tal emission is seen up to 2" (240 A U projected distance) from 
the central source. In a larger sample, Pontoppid an et alJ (2003) 
detected bright CO emission from 18 of 44 objects observed with 
VLT/ISAAC at moderate resolution (R = A/ 5 A = 10 4 ), but did 
not analyze this emission because, when spectrally unresolved, 



multiple absorption and emission components can introduce 
spurious velocity shifts and cause misleading non - detec tions 
(see also iPontoppidan et all l2005al) . IBrittain et all d2005l) at- 
tributed CO emission from the embedded object HL Tau to the 
disk, based on similarity of the emission profile seen at high res- 
olution to the profiles from CT TSs that lack envelopes . Early 
studies bv lChandler et all d 19931) and iNajita etafl d 19961) found 
vibrationally-excited CO overtone (v = 2 - 0) emission from two 
embedded objects, WL 16 and LI 551 IRS 5. 

Now that CO v = 1 - emission has been used to study large 
samples of disks around CTTSs and Herbig AeBe stars, apply- 
ing the same technique to younger YSOs could reveal the pres- 
ence and structure of young disks that are difficult to detect with 
other gas tracers. Embedded objects should have larger accretion 
rates than classical T Tauri stars (CTTSs), which may introduce 
differences between young disks around embedded objects and 
more mature disks around CTTSs. Any warm (500 - 2000 K) 
molecular gas in the inner disks of Stage I YSOs should produce 
bright and broad CO emission lines. CO fundamental emission 
from warm gas has the potential to be a clean diagnostic of disks 
during the embedded phase of YSO evolution and is free from 
contamination by colder emission expected from both the outer 
disk and envelope. Moreover, even if a compact, disk-like struc- 
ture is present, as inferred from sub-mm interferometry, resolved 
line profiles could indicate w hether the gas has h ad time to settle 
into Keplerian rotation (e.g. lBrinch et a l.. 2009). 

In this paper, we explore the utility of CO fundamental emis- 
sion as a possible probe of these young disks by analyzing high 
spectral resolution (R = A/AA ~ 90,000) M-band observations 
of 1 8 YSOs observed with CRIRES on the Very Large Telescope. 
CRIRES is uniquely capable of high-resolution M-band spec- 
troscopy with relatively broad spectral coverage at high spatial 
resolution. The data presented here were obtained as part of 
a large program to survey CO emission fr om YSOs at d iffer- 
ent sta ges of pre-main sequence evolution (Pontop pidan et all 
1201 lbl) . Most observations in this paper were obtained under ex- 
ceptional weather conditions. In §2, we describe the observa- 
tions and the selected sample. In §3, we describe and analyze 
the observed line profiles, including separating the emission into 
broad and narrow components, analyzing spatially-extended CO 
and H2 emission, and finding CO absorption in winds of six ob- 
jects. In §4, we attribute the components to different regions in 
the YSO and discuss the implications of our results on disk struc- 
ture and wind launch regions. 



2. OBSERVATIONS 

2.1. Target Selection 

The selected targets are Stage I 1 YSOs that are still embedded 
in circumstellar envelopes and that are members of nearby star- 
forming regions. Most of the targets in our sample were cho- 
sen based on fundamental CO emissio n that was previously de- 
tected in R ~ 10, 000 ISAAC spectra dPontoppidan et all 120031 
l2005ah . Most targets are low-mass YSOs based on their lumi- 
nosity. One object, SVS 20 S, may have a luminosity consistent 
with intermediate-mass YSOs. A few other objects may evolve 
into Herbig AeBe stars. 

Because SEDs of heavily-extincted stars can be similar to 
SEDs of stars embedded in circumstellar envelopes, the presence 



1 For simplicity, we use the physical terminology "Stage" rather than 
the observational terminology "Class" throughout the paper, although 
in some contexts "Class" would be the more appropriate term. 
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Table 1. Evidence for envelopes within our sample 



Source 


mm-int." 


Compact 6 


Outflow'' 


Ref 


TMC 1A 


Y 


Y 


Y 


1,2,3,7 


GSS 30 


Y 


Y 


n 


1,2,4 


WL 12 


Y 


Y 


s 


1,2,4 


Elias 29 


Y 


Y 


Y 


1,2,4,8 


IRS 43 S+N 


Y 


Y 


s 


1,2,4 


IRS 44 E+W 


- 


Y 


Y 


1,4 


IRS 63 


Y 


Y 


Y 


1,2,4 


L1551 IRS 5 AB e 


Y 


Y 


Y 


3,7,13 


HH 100 IRS 




Y 


Y 


5,9 


WL 6 




Y 


n 


1.4 


Elias 32 




Y 




4 


CrA IRS 2 










HLTau 


Y 


Y 


Y 


3,10,12 


SVS 20 S+N 




Y 




7 


RN0 91 




Y 


Y 


6,11 



"Envelope detected in mm-interferometry. 

^Presence of cold, compact gas centered at object position. 

^Refers to cold molecular outflows detected in the (sub)-mm. 

(s): Molecular outflow is detected in only a single lobe. 

(n): No molecular outflow detected. 

''The references are examples of each phenonemon and 

are not intended to be complete. 
''Unresolved and hereafter referred to as L1551 IRS 5 
1: Bontemps et al. (1996) . 2: Jotgensen et al. (20091 
3 : Saito et al. (2001) . 4: van Kempen et al. (2009) 
5: van Kempen et al. (2009). 6: Arce et al. (2006) 
7: Gregersen et al. (2000) . 8: Lommen et al. (2008) 
9: Groppi et al. (2007) . 10: Wilneretal. (1996) 
11: Chenetal. (2009) . 12: Cabrit et al. (1996) 
13: Wu et al. (2009) 



of an envelope for objects in our sample is typically confirmed 
by the presence of at least one of the following three criteria: (1) 
a clearly exte nded component in the visibility curves of (s ub)- 
mm emission ( Lommen et al .1 12008t I.Tdmensen et all l20Q9h . (2) 
a compact structure of c old dense ga s , such as HCO + , at the same 
position as the targe t (ISaito et all 120011: iGroppi et all l2007t 
lLommen et ail l2008t Ivan Kempen et all \2QQ% . or (3) a cold 
molecular outflow se en close to the YSO in both a red- an d 
blue-shifted lobe (e.g. lBontemps et all [19961 lArce et aill2006h . 
Table Q] summarizes the evidence for our classification. One ob- 
ject, CrA IRS 2, is less-well studied and not confirmed as a Stage 
1 object by these cr iteria, but is likely embedded in an envelope 
dNutter et all 120051) . For the purposes of this paper, both ob- 
jects in spatially-resolved multiple systems with envelopes are 
assumed to be embedded (see Appendix A for more details). 

The classification of these targets as embedded is consistent 
with most previous classifications, but can differ from classifica- 
tions t hat are based on SEDs alone. For example. iMcClure et alj 
(l2010h used 12-5 yum colors to classify WL 12, IRS 43, and IRS 
44 as YSOs with envelopes and WL 6, Elias 29, and IRS 63 
as YSOs with disks but no significan t envelope. When c lassi- 
fying the latter three objects as disks, McCl ure et all (1201 Ol) ex- 
plain the discrepancy with previous work by suggesting that the 
envelopes are tenuous and close to d isappearing. This interpre - 
tion may apply to WL 6 and IRS 63 (van Kemp en et all 12009). 
However, the sub-mm visibility of Elias 29 shows no indica- 
tion of any compact structure, which means that the envelope 
mass must be much la rger than the mass of the undetected disk 
(lLommen et all 120081) . The disks for these three sources may 
dominate their near-IR spectra, but the sub-mm interferometry 
requires the presence of an envelope. 



Table [2] lists selected properties of our targets from the 
literature. The stellar properties, including the effective tem- 
perature and luminosity of the photospheric emission and 
the extinction, are uncertain bec ause of the difficulty in de- 
tecting photospheric fe atures (e.g. IWhite & Hillenbrand [2004; 
Doppman n et all 1 20 05) and in accurately measuring the extinc- 
tion. In some cases, the Ay and Lb \ are obtained from different 
sources. The bolometric luminosities that are measured from the 
far-IR/sub-mm SEDs are not particularly sensitive to either ex- 
tinction or spectral type. All luminosities are corrected for the 
most recent distance measurement to the parent cloud. In a few 
cases, significant discrepancies in derived stellar parameters ex- 
ist in the literature 2 . 

The sample likely includes a range of masses for the central 
object, with a general trend that more luminous objects are likely 
also more massive. Unfortunately, the lack of accurate stellar 
properties prevents a reliable determination of mass from pre- 
main sequence evolutionary tracks. 

The systemic velocity listed in Table|2]is obtained from liter- 
ature measurements of (sub)-mm lines of the molecular cloud or 
of the envelope. All CO and H2 velocities listed in the paper are 
relative to the measured velocity of I3 CO and C I8 absorption 
in our M-band spectra. 

2.2. Observational Setup and Data Reduction 

We used CRIRES (iKaufl etalll2008h on the VLT-UT1 to obtain 
high-resolution echelle spectra from 4.6-4.9 //m of 15 Stage I 
objects, three of which are binaries. Table[3]lists our observation 
log 3 . 

CRIRES has four 1024 x 512 pixel InSb detectors that each 
covers ~ 0.02 ^m in a single integration. Every object in our 
sample was observed at multiple wavelength settings to cover 
spectral gaps between chips and to observe lines with a wide 
range of rotational levels. Table [3] includes the wavelength set- 
ting used for each object. Online Table A. 1 lists the wavelength 
ranges covered for each setting. Each pixel covers ~ 0.2 A, or ~ 
1.25 km s _1 , in the dispersion direction. The 0719-wide slit with 
07086 pixels in the dispersion direction leads to R ~ 90, 000, as 
measured from telluric absorption lines. 

Observations were obtained with 10" nods in an ABBA pat- 
tern. Each nod consists of 60s integrations. Random dither pat- 
terns with a width of 1" were used for each nod to distribute the 
counts over different pixels at each integration. Total integration 
times for each setting are listed in Table [3] 

Most sources in our sample were too faint for the MACAO 
adaptive optics system, which usually feeds CRIRES. The AO 
loop was closed for HL Tau and RNO 91, though with a low 
Strehl ratio and poor correction. For other objects, the M-band 
seeing during our observations was exceptional, typically be- 
tween 073 and 0745 as measured from the spatial profile of the 
M-band continuum emission in the cross-dispersion direction for 

2 For IRS 63, th e photospheric properties measured by 
Doppmann et al. ( 2005|) are not adopted here because the v; sr = -26.5 
km s"' of photospheric absorption features is highly anomalous 
relative to the other Stage I objects in the p Oph molecular cloud and 
is inconsistent with the velocity o f CO absorption by th e envelope 
seen in our data. We speculate that Dopp mann et alj d2005l) measured 
optically-thin absorption in an outflow and attributed the absorption to 
the photosphere. Alternately, the slit rotation may have induced a large 
velocity offset, in which case the derived effective temperature would 
still be accurate (T. Greene, private communication). 

3 Data are available for download at 
http://http://www.stsci.edu/~pontoppi 
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Table 2. Sample Properties 



Target 


d 


PA" 


A v 


T piwt 


Lboi 


v; sr 


abs 


Flux'' 






pc 


° 


mag 


K 


Lo 


km s" 1 


km s' 1 


4.8 yum 


Method 


Ref 


HH 100 IRS 


130 


- 


30 


4060 


15 


5.9 


5.8 


11.3 


ISO 


2,8 


CrA IRS 2 


130 


- 


20 


4900 


12 


(5.9) 


6.1 


5.40 


ISO 


8,15 


WL6 


120 


-40 


9.8 


- 


2.4 


3.5 


4.5 


1.4 


IRAC 


1,23 


WL 12 


120 


- 


9.8 


4000 


2.4 


3.5 


7.6 


1.1 


IRAC 


1, 10 


IRS 63 


120 


~ 240 


9.8 


- 


3.0 


3.5 


2.0 


1.3 


IRAC 


1,24 


IRS 43 S+N 


120 


25 


9.8 


4400 


5.5 


3.5 


4.1 


1.6 


IRAC 


1,10,17 


IRS 44 E+W 


120 


-0 


9.8 


4400 


25 


3.5 


5.8 


2.8 


IRAC 


1,10,20 


Elias 32 


120 


- 


9.8 


3400 


1.0 


3.5 


5.7 


0.39 


IRAC 


1,10 


Elias 29 


120 


~ 135 


9.8 




38 


3.5 


5.5 


24.3 


ISO 


1,21 


GSS 30 


120 


-45 


9.8 




13 


3.5 


7.5 


23 


IRAC 


1,20 


HLTau 


160 


-45 


7.4 


4350 


6.6 


7 


8.2 


5.49 


ISO 


3,4,9,22 


RN0 91 


120 


155 






2.3 


0.8 


0.5 


1.8 


16 


5,12,18 


SVS20 S 


415 




30 


5900 


142 


1.9 


8.3 


3.5 


IRAC 


10,14,7 


SVS20 N 


415 




4 


3270 


0.27 


1.9 


(8.3) 


1.1 


IRAC 


13,7 


TMC 1A 


160 


-10 






2.8 


5.6 


5.3 


1.0 


IRAC 


6, 3,7, 25 


L1551 IRS 5 AB 


160 


256 


28 


3300 


23 


6.9 


4.7 


3.94 


ISO 


6, 3,7,10,11,19 



"Approximate position angle of outflow *v/„ of 1J CO and C 1!i O absorption in our CRIRES spectra 

c 10" 14 erg citT 2 s A -1 () indicates assumed value from nearby source. 

Distances: 130 pc to CrA (de Zeeuw et al.. 1999) . 120 pc to Oph (Loinard et al.. 2008) . 415 pc to Serpens 



(Dzib et al.. 2010) . and 160 pc to Taurus (Tor res et al.. 2009) 
1: Evans et al. (2009) . 2: van Kempen et al. (2009) 3 :Furlan et al. (2008) 4: Havashi et al. (1993) 
5: Lohretal. (2007) : 6: Yang et al. (2002) 7: Gregersen et al. (2000) 8: (Nisini et al.. 2005) 
9: White & H illenbrand (2004) : 10: Doppmann et al. (2005) : 11: Pratoetal. (2009) : 12: Chen et al. (2009) 
13: Oliveira et al. (2009) : 14: Ciardi et al. (2005) : 15: Mever et al. (2009) : 16: Booge rt et al. (2008) 
17: Grosso et al. (2001) : 18: Arce et al. (2006) : 19: Pvo et al. (2009) : 20: Allen et al. (2002) : 21: Ybarraetal. (2006) 
22: Beckwith et al. (1989) : 23: Gomez et al. (2003) : 24: Zhang & Wang (2009) : 25: Chandler et al. (1996) 



Table 3. Observation Log 



Target 


Alt. Names 


RA" 


DEC" 


UT Date 


tap (S) 


airmass 


FWHM" 


Vbary (km s ') 


A settings 


S/N" 


IRS 43 


YLW 15 


16:27:27 


-24:40:51 


2008-08-06 


960' 


1.06 


0.32 


-17 


4716, 4730, 4868 


40 


IRS 43 


YLW 15 


16:27:27 


-24:40:51 


2008-08-07 


1200 


1.02 


0.25 


-17 


4831 


40 


IRS 44 rf 


YLW 16a 


16:27:28 


-24:39:34 


2008-04-27 


720 


1.41 


0.56 


27 


4716 


15 


IRS 44 


YLW 16a 


16:27:28 


-24:39:34 


2008-04-30 


720™ 


1.26 


0.54 


25 


4730, 4833 


10 


IRS 44'' 


YLW 16a 


16:27:28 


-24:39:34 


2008-05-01 


1200 


1.61 


0.86 


25 


4946 


9 


IRS 44" 


YLW 16a 


16:27:28 


-24:39:34 


2008-08-06 


1200 


1.31 


0.33 


-17 


4868 


20 


IRS 44'' 


YLW 16a 


16:27:28 


-24:39:34 


2008-08-07 


960 


1.38 


0.32 


-17 


4716,4730 


30 


IRS 63 


GWAYL 4 


16:31:36 


-24:01:29 


2007-04-26 


600 


1.02 


0.37 


28 


4716, 4730, 4833,4929 


32 


WL 12 


YLW 2 


16:26:44 


-24:34:48 


2007-09-02 


960 


1.22 


0.31 


-20 


4716,4730 


48 


WL 12 


YLW 2 


16:26:44 


-24:34:48 


2007-09-04 


960 


1.25 


0.31 


-20 


4833 


45 


WL 12 


YLW 2 


16:26:44 


-24:34:48 


2010-03-03 


2280 


1 .05 


0.31 


40 


4946 


60 


WL 6 


YLW 14 


16:27:22 


-24:29:53 


2007-04-27 


480 


1.20 


0.29 


27 


4716, 4730, 4833 


30 


HH 100 IRS 


V710 CrA 


19:01:51 


-36:58:10 


2007-09-03 


600 


1.27 


0.43 


-17 
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"2MASS position (12000) from (Skrutskie et al.. 2006). 

"In arcsec, includes both the seeing and any real spatial extension in the continuum. 

"Approximate S/N per pixel in a single setting, typically measured in regions with an average telluric absorption of 10-15%. 
Position angle of 0° except where noted: d PA=70; «PA=100; ^PA=34(); «PA=-150; ''PA=45; 'PA=5; JPA=160; *PA=75; 
Listed position angle is approximate for eases where the slit was aligned with the parallactic angle. 

Exposure times same for all settings except where noted: '4868: 1440s; "'4833: 960s "4868: 960s "4868: 4929/4946: 500; "4946: 1080; ''51 15: 1080s 
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Fig. 1. Extraction of two components from the IRS 44 spectra. 
The main plot shows the extracted counts versus wavelength for 
IRS 44 E (red) and IRS 44 W (blue) over a small wavelength 
range. The two insets show the spatial profile of the counts in 
the cross-dispersion direction for a single pixel (left, within the 
CO P(21) line, right within the continuum). Two Gaussian pro- 
files were fit to the cross-dispersed emission profile and were 
subtracted to minimize contamination between components dur- 
ing the extraction. CO emission is detected from the secondary, 
IRS 44 W (blue) but not from the primary, IRS 44 E (red). The 
excised region is unreliable because of strong telluric CO ab- 
sorption. 



each object. Two objects, Elias 29 and L1551 IRS 5, have espe- 
cially large seeing measurements, which indicates that either the 
seeing was poor or that the M-band continuum emission is spa- 
tially extended on 0/5 - 170 scales. 

Each object was acquired with the Ks filter. After November 
2008, the CRIRES control software implemented an automatic 
correction to the differential atmospheric dispersion between the 
K and M-bands. Prior to November 2008, our observations were 
generally obtained at low airmasses or at parallactic angle to 
minimize the effect of atmospheric dispersion. Because most of 
our observations did not use adaptive optics, the effect of atmo- 
spheric dispersion (~ 50 mas at airmass of 1.7) is typically much 
less than the FWHM of the continuum in the cross-dispersion 
direction. A real offset between the K-band and M-band con- 
tinuum emission is also possible for embedded young stars. In 
cases where the slit was aligned with the parallactic angle, the 
position angle did not change by more than 10° in any observa- 
tion. All of our observations include strong continuum emission, 
and effects that may be attributed to different slit placements are 
likely minor. 

A telluric standard at a similar airmass to our object spec- 
trum was always observed in an adjacent observation and with 
an AO correction. Spectra from the object and the telluric stan- 
dard were wavelength-calibrated using telluric absorption lines. 
The science spectrum was then divided by the spectrum of the 
telluric standard. Differences in spectral resolution between the 
telluric standard and the science target were always minor and, 
when necessary, were accounted for by convolving the spectrum 
of the telluric standard with a Gaussian profile to match the two 
resolutions. Figure [2] demonstrates that high S/N (~ 50 - 100 
in the continuum) is obtained even in regions filled with telluric 
absorption lines. The spectra were then shifted in velocity space 
to the local standard of rest. Many of our observations were ob- 




4.650 4.655 4.660 
Wavelength (fim) 



4.665 



Fig. 2. The telluric correction for IRS 43 S+N from two chips in 
the 4716 A setting. The spectrum reaches S/N ~ 100 per pixel 
in regions with few atmospheric lines (left) and S/N ~ 50 per 
pixel in regions with many atmospheric lines (right). I2 CO lines 
are labeled in black and 13 CO lines are labeled in red. 



tained on dates when the telluric absorption was shifted by large 
velocities from the local standard of rest to shift telluric CO ab- 
sorption away from the central wavelength of CO lines for the 
objects. The relative wavelength calibration is accurate to ~ 1 
km s _1 . 

Data were reduced following lPontoppidan et all d2008l) . The 
pixels on the detector are each O/086 in the dispersion and cross- 
dispersion direction. When coadding the 2D images, the light 
was shifted onto a common central pixel and resampled onto a 
finer (07043) scale. Spectra were extracted from windows with 
a width that is twice the FWHM of the emission profile in the 
cross-dispersion direction. A specialized extraction technique 
was applied to the 0/3 binary IRS 44 (see §A1.1), which was 
marginally resolved on three nights when the seeing was excep- 
tional (073). The spectra for the two components were extracted 
by fitting two Gaussian profiles to the cross-dispersed profile at 
each spectral bin in the 2D images (see Fig . [TJ . The two Gaussian 
profiles both have the same FWHM. The FWHM and positional 
offset between the two components are constant across the detec- 
tor. The Gaussian profile of one component was subtracted from 
the data. The counts were then extracted from a spectral win- 
dow around the other component. The extraction windows for 
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Fig. 3. A sample spectral region, with 12 CO v = 1 - lines shown in black, 13 CO v = 1 — shown in red, C ls O v = 1 — shown in 
blue, and 12 CO v = 2 — 1 shown in green. 



the two components do not overlap to minimize contamination 
by the other component. 

The M-band continuum flux i s measured d irectly from ISO- 
SWS spectra for four objects (Sloan et al 1 12001 and from 
groun d-based ISAAC spectra for one object (Boogert et al., 



2008), and is inferred f rom linear fits to SpitzerfLRAC photom- 
etry f or seven objects (Ivan Kempen et al.L 120091; lHarvev et all 
2007). For binaries, the continuum flux of each companion was 
calculated by using the M-band magnitude and the counts ra- 
tio in the continuum emission measured in our M-band observa- 
tions. Spatially-extended continuum emission for a given object 
would artificially increase the calculated M-band flux of the cen- 
tral source. The CO line fluxes are calculated from equivalent 
width measurements and would be overestimated in cases where 
the M-band continuum is spatially extended. 

Our CRIRES spectra yield a flat M-band continuum flux for 
every object, with a relative uncertainty of ~ 30% between 4.5 
and 5.0 fim. In the four ISO spectra, the maximum relative flux 
difference between 4.5 and 5.0 /mi is 30%. We also do not cor- 
rect the absolute or relative CO emission line fluxes for extinc- 
tion because of the wide range in calculated extinctions for any 
given object, the lack of methodological consistency in extinc- 
tion calculations for the entire sample, and because extended CO 
emission may not suffer from the same extinction as the cen- 
tral source. Applying extinctions that range from Ay = 10 - 37 
mag would increase CO l i ne luminosities by fac tors of 1.2 to 3 
dRieke & LebofskvL 119851: IChapman et ail [2009) relative to the 
measured luminosities reported in this paper. The wavelength 
dependence of extinction in the M-band is negligible (~ 3% 
across the 4.6-5.0 fj.m spectral region for Ay = 20 mag). 



2.3. Archival ISAAC Spectra 

In several subsections we compare our high-resolution CRIRES 
spectra to archival VLT/ISAAC M-band spectra of the same 
sources obtained in 2001-2002. The ISAAC spectra span from 



4.53^1.75 yum with R = 5000 - 10000. The spectra were ob- 
tained with a 073 slit width and generally with worse seeing 
than our CRIRES data, possibly because the embedded ob- 
jects within the CRIRES program were preferentially observed 
on nights with good seeing. As a consequence, the extracted 



ISAAC spectra usually sample lar ger areas on the sky. The data 
were reduced and analy zed by IPontoppid an et al.l d2002l) and 
iPontoppidan et alj (2003), including reporting detections and ve- 



locities for CO emission. 



2.4. Archival Near-IR Images of GSS 30 

To compare extended CO emission to continuum emission from 
GSS 30 (see §3.3), we downloaded an AO-fed K-band image 
from the VLT archive that was obtained with NACO on 19 
September 2007 as part of VLT program 079.C-0502 (P.I. Chen). 
Gaspard Duchene (private communication) provided us with an 
L-band AO image of GS S 30 o btained with VLT/N ACO and pub- 
lished in lDuchene etall (l2007l) . 



3. DESCRIPTION OF M-BAND EMISSION FROM 
EMBEDDED OBJECTS 

M-band spectra cover many CO fundamental (Av =1) transi- 
tions, H Pfund-/?, H 2 S(9), and the 4.67 /im CO ice absorption 
band. One spectral region of our CRIRES M-band spectra of 
embedded objects is shown in Figure [3] The M-band emission 
is dominated by a smooth continuum produced b y warm dust, 
likely located in the dis k close to the star (e.g., lEisner et"aU 
120051 : lEnoch et all 12009b . No photospheric f eatures are detected 
from any of our targets. CO ice absorption (Pontop pidan et all 
2003) is detected towards most objects in our sample. HL 
Tau lacks any CO ic e absorption (see also IWhittet et all fl989: 
Brittai n et all 120051) . which suggests that much of the line- 
of-sight extinction occurs in a disk or flattened envelope in 
which little solid CO absorption is expected to be observed 



6 



Herczeg et al.: CO Emission from Embedded Objects 




-150 -75 75 150 

Velocity (km s -1 ) 



Fig. 4. Co-added profiles for CO v = 1 - 0, v = 2 - 1 and 13 CO 
v = 1 - lines in the VL77CRIRES spectrum of HH 100 IRS. 
The different observed components are each identified with a 
physical region in the system. For scaling, the peak of the v = 
1-0 line is at ~ 1.38, normalized to the adjacent continuum. 



dPontoppidan et all l2005bh . Strangel y, GSS 30 also shows only 
very weak ice absorption (see also Pontopmdan et al., 2002). 
Absorption in gaseous 12 CO and isotopologues is detected to- 
wards every object. Narrow absorption lines are caused by the 
envelope and parent molecular cloud, while broad absorption 
arises in the wind. 

Emission in 12 CO v = 1 — lines is detected from all but 
four embedded objects in our sample. In general, the on-source 
CO emission can be described by a broad and a narrow com- 
ponent. Broad CO emission, in both v = 1-0 and in higher 
vibrational lines, is detected from 10 sources in our sample. The 
narrow component, which has an optical depth sufficient to pro- 
duce detectable emission in 13 CO and C 18 lines, is detected 
from 9 sources. These different components are more easily dis- 
tinguished in the v = 2 - 1 and isotopologue lines that are also 
covered in M-band spectra. Figure [4] shows how the broad com- 
ponent, seen in v = 2 - 1 lines, and the narrow component, seen 
in 13 CO lines, combine to create the 12 CO v = 1-0 emission 
line profiles for one object, HH 100 IRS. Similarly for the other 
sources, the I3 CO lines are typically narrower than the v = 2 - 1 
lines (Fig. [3]). In addition to the on-source emission, very nar- 
row CO emission is spatially-extended from two sources (GSS 
30 and IRS 43, see §3.3). Blueshifted CO absorption is detected 
from six sources (see §3.6). 

Figure[6]presents a gallery of co-added CO line profiles from 
all objects with detected CO emission. Lines were co-added over 
all J transitions that were observed and that have minimal or 
moderate contamination from telluric absorption. Some objects 
show strong emission in the narrow, optically-thick component, 
some objects show strong emission in broad, vibrationally ex- 
cited emission, other objects show both, and still other objects 
show only weak or no emission in 12 CO v = 1 — lines. Figure|7] 
and Table [4] (ordered by decreasing luminosity) summarize the 
presence and absence of each component for these objects. 
Objects with a high bolometric luminosity, relative to the median 
in our sample, show emission in the narrow component but not 
the broad component. Objects with a lower bolometric luminos- 
ity, relative to the median in our sample, emit in the broad com- 
ponent but not in the narrow component. Hereafter, we loosely 
define high and low bolometric luminosities as greater than or 
less than ~ 15 L G , respectively. 



Table [5] summarizes the properties of Gaussian profiles that 
were fit to the different isotopic and vibrational lines. The fits 
were applied to lines coadded over all unblended rotational lev- 
els (Fig. [6]). The emission from v' > 1 is typically broad (FWHM 
of 50-160 km s _1 ) and centered at the systemic velocity. None 
of the lines show a double-peaked profile that is characteristic of 
Keplerian rotation in a disk with high inclination, although the 
co-added CO v = 2 — 1 profiles from RNO 91 tentatively show a 
weak dip in flux relative to a Gaussian profile. On the other hand, 
the 13 CO emission lines are narrow (FWHM between 10 to 50 
km s ), with central velocities shifted between to -15 km s . 
The range of velocity centroids for the narrow component sug- 
gests that our crude classification scheme does not capture the 
full range of physical processes that produce this component. 

In most cases, the CO emission can be clearly distinguised as 
a broad or narrow component. In a few spectra, the precise clas- 
sification of CO components is unclear. For example, the narrow 
component of CO emission from GSS 30 and IRS 44W each in- 
clude emission from at least two regions based on distinct com- 
ponents in the line profiles. These cases are discussed in more 
detail in the following subsections and in Appendix C. 

Online Table A. 2 presents the line fluxes and equivalent 
widths for the different components in several selected transi- 
tions. Most 12 CO v = 1 -0 and 13 CO v = 1 -0 lines are partially- 
attenuated by absorption in the same transition. The listed line 
fluxes are what the fluxes would be if no line-of-sight CO ab- 
sorption occurred. The fluxes are not corrected for extinction. 
The absorbing CO gas in the circumstel lar envelope or molec- 
ular cloud is often cold (~ 10 - 100 K. |j0rgensen et all 120021: 
Bergin & TafaliaL l2007l) and only affects emission in low-/ tran- 
sitions. However, some absorption components are optically- 
thick for transitions with / > 30. In some cases, multiple CO 
absorption components are detected. The listed fluxes are mea- 
sured from Gaussian fits to the emission profile, ignoring any 
absorption component. The FWHM and central velocity are de- 
termined from fits to the summed profiles (Tabled. The central 
velocities are listed relative to the velocity of 13 CO and C 18 ab- 
sorption listed in Table [2] Fitting profiles to the observed emis- 
sion while ignoring the absorption components yields line fluxes 
measured consistently throughout a given spectrum. In many 
cases the fit is applied only to a portion of the line profile to 
avoid absorption features, regions with optically-thick telluric 
absorption, contamination from emission in other components 
of the line, and emission in other lines. In each case the 1-cr 
uncertainty in flux is dominated by the uncertainty in the contin- 
uum level and does not include the uncertainty in the relative or 
absolute flux calibration. 

In the subsequent subsections, we describe separately the 
properties of the broad and narrow emission components, 
spatially-extended CO emission, H2 emission, and CO wind ab- 
sorption detected within our sample. For each CO emission com- 
ponent, a model of an isothermal, plane-parallel slab of CO gas 
(see Appendix D) is used to calculate temperatures, column den- 
sities, and emitting areas, with results in Table[6] 



3.1. Broad CO Emission From Warm Gas 

Broad emission in v = 2 - 1 lines is clearly detected in 10 of 18 
sources 4 (Figure |5}. Emission in some v = 3 - 2 and 4-3 lines 
is detected from CrA IRS 2 and HH 100 IRS, which have high 
S/N. Non-detections of lines from V > 2 from other sources are 
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Fig. 5. Line profiles of 13 CO (left) and CO v = 2 — 1 (right) emission for detected sources in our sample. For comparisons to CTTSs, 
the v = 2 — 1 line profiles of Elias 23 and S CrA S and the 13 CO line profiles of S CrA S and DR Tau are shown at the top. 



Table 4. Presence of CO emission components 
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generally not significant, assuming the same flux ratio in v 
3-2tov = l- lines as in CrA IRS 2 and HH 100 IRS. 



4 Counting WL 12 as a detection, despite a non-detection in one of 
the two observations, and GSS 30 as a non-detection, despite some v = 
2-1 emission (see §3.2.1). 



Fig. 7. A Venn diagram showing the presence of CO emission 
and absorption components from Table [4] 



In most cases, a component of similar shape and strength as 
the v = 2 — 1 lines is seen in the v = 1 - lines. The v = 
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Fig. 6. Co-added line profiles from sources with detected CO emission. In each plot the 12 CO v = 1 — line is shown in red, the 
13 CO v = 1 - line is shown in blue, the C ls O v = 1 - line is shown in pink (IRS 44 W and GSS 30), and 12 CO v = 2 - 1 line 
is shown in green. Inserts zoom in to compare the 12 CO emission line profile to that of 13 CO, C ls O, or v = 2 - 1, depending on the 
source. 



1-0 emission is often absorbed by CO in our line of sight to the 
emission region. Broad CO wind absorption partially or totally 
obscures any CO v = 1 — emission on the blue side of the 
line profile of six sources (HH 100 IRS, IRS 63, Elias 29, TMC 



1A, Elias 32, and WL 12, see §3.7). In most cases, the same 
broad component seen in v = 2 — 1 lines is also identified in 
the v = 1 - lines based on the similarity of the emission line 
profiles on the red wing. For HL Tau and CrA IRS 2, CO v = 
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Table 5. Average CO Emission Line Properties 
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18(7) 


Elias 29 e 


v = 1 - 


(3.0) d 


19(3) 


IRS 63" 


v = 1-0 


3(10) 


28 (13) 


HLTau 


13 co 


-15 (6) 


42(14) 


SVS 20 N 


v = 1 - 


-8.8 (1) 


28 (5) 


CrA IRS 2 f 


13 co 


0.0 (0.4) 


34 (4) 


CrA IRS 2 f 


v = 1-0 


0.5 (0.3) 


16(3) 


SVS 20 S 


v = 1-0 


-8.0(1.0) 


37 (7) 


SVS 20 S 


13 co 


(-8.0) rf 


47 (23) 


IRS 44 W 


13 co 


15.4 (0.4) 


14(1) 


IRS 44 W 


c 18 o 


(15.4) d 


8(2) 


IRS 44 W 


12 co 


-1.0(1.0) 


32 (4) 


L1551 IRS 5 


13 co 


0.8 (0.7) 


12(2) 


L1551 IRS 5 


12 co 


(0.8) rf 


17(2) 


L1551 IRS 5 


v = 2 - 1 


(0.8y 


13(2) 


GSS 30'' 


v = 2-l 


-5.8 (2.0) 


42 (10) 


GSS 30 


v = 2-l 


-19.2 (0.5) 


10(2) 


GSS 30 /! 


12 CO,/ > 30 


-10(1) 


36 (3) 


GSS 30 SW 


12 CO 


-9.1 (1.0) 


6.3 (1.0) 


GSS 30 NE 


12 CO 


-8.6(1.0) 


6.5 (1.0) 


IRS 43 S'' 


v= 1 -0 


-4(1) 


42 (7) 



CO transition 



"Based on Gaussian profiles fit to coadded lines. 
^Relative to velocity of l3 CO and C 18 absorption. 
Trom spectrum obtained in August 2008 
"Forced value from same component in v = 2 - 1 or 13 
Tit only to red emission 
Tit is not good because line is not Gaussian 
s March 2010 

'Tentative classification as narrow component 



Table 6. Physical Properties of CO Emission Region 



Star 


Component 


Tro, (K) 


logW(CO)" 


(Area) " , 


IRS 44 W 


Narrow'' 


330±20 


19.15±0.25 


3.6 , 


GSS 30 


Narrow'" 


315 ±15 


18.65 ±0.10 


fl ' 


GSS 30 


Extended 


250 ± 30 


18.7 ±0.2 




CrA IRS 2 


Broad 


1100 ±200 






CrA IRS 2 


Narrow 


560 ± 70 


19.1 ±0.1 


1.0 ' 


HH 100 IRS 


Broad 


1000 ± 100 


17.8 ±0.2 


0.7 ; 



"Column density N in units of log cm - - here 

and throughout paper, assumes b = 2.0 km s - 
''Square root of emitting area, in units of AU 
Tlueshifted emission 



1-0 emission in a narrower emission component masks any 
possible broad emission that would have the same profile as the 
v = 2 - 1 emission. The broad emission component is not seen 
in any 13 CO lines within our sample. 

The centroid of the v = 2 - 1 emission is usually consistent 
with the systemic velocity. The centroids that deviate by > 3 
km s _1 (IRS 43 S, HL Tau, and WL 12) each have very low S/N 
in the coadded line profiles and as a result have unreliable central 
velocities. In the cross-dispersion direction, the broad emission 
from HH 100 IRS, RNO 91, and CrA IRS 2 (the three cases with 
highest S/N in the broad lines and the highest spatial resolution) 
is centered at the same location (to within ~ 0.2 pix, or 1-2 AU at 
120 pc) as the continuum emission and is not spatially extended 
(FWHM< 1 - 1 .5 pix, or 6-8 AU at 120 pc). 

To measure temperatures and column densities (Fig. [8), we 
select the two spectra, HH 100 IRS and CrA IRS 2, that have the 
highest S/N in broad line emission over a wide range of / levels. 
The excitation diagram of CO emission from HH 100 IRS shows 
an upturn in N/(2J+ 1) at low-/, characteristic of optically-thick 
emission from warm (500-1500 K) gas (see Appendix D). The 
CO v = 2 — 1 lines are more likely to be free of any optical 
depth effects. No 13 CO v = 1 — emission is detected, with an 
upper limit that is ~ 10 times weaker than the 12 CO v = 1 - 
emission. A fit to the v = 2 - 1 lines indicates a rotational tem- 
perature of ~ 1000 + 100 K. At this temperature, the lack of 
detectable 13 CO emission in this component requires a column 
density logA^(CO)< 18.0. In the excitation diagram, the curve 
for v = 1 — lines requires logA^(CO)> 17.5. The emission 
in 12 CO v = 1 - lines with high-/ is underproduced at this 
temperature and column density, which may suggest that the ro- 
tational temperature derived from the v = 2 — 1 lines is too low 
or that a second temperature is needed. The v = 2 — 1 lines are 
underproduced by ~ 0.6 dex, indicating that the vibrational tem- 
perature is warmer than the rotational temperature. The flux ratio 
of v = 3 - 2 and 2-1 lines yields a vibrational temperature of 
~ 2000 K. The total emitting area is roughly (0.7 AU) 2 . 

For CrA IRS 2, the broad component in the v = 1 — lines is 
not measureable because it is mostly masked by much stronger 
emission in a narrower emission component (see Appendix C.3). 
Broad emission is detected in v = 2 - 1, 3 - 2, and a few 4-3 
emission lines. The v = 2 — 1 and 3-2 line flux ratios yield a 
rotational temperature of ~ 1000 K, while the ratio of line fluxes 
in the different vibrational levels yields a vibrational temperature 
of ~ 2500 K. 



3.2. Narrow CO Emission From Optically-Thick Gas 

Emission in 13 CO v = 1-0 lines is detected from 9 of 18 
embedded objects within our sample. The same narrow compo- 
nent seen in 13 CO emission is also seen in 12 CO and C 18 lines, 
which indicates that the 12 CO lines are optically-thick. Some of 
.these same objects also have a broad emission component. In 
only one case, GSS 30, is CO v = 2 — 1 emission detected in 
the narrow component, probably because of the high S/N and 
line-to-continuum contrast in that spectrum. 

The line profile of narrow emission is centered at the sys- 
temic velocity for most sources. In addition to emission at the 
systemic velocity, IRS 44 and GSS 30 also show narrow com- 
ponents blueshifted by ~ 10 km s _1 . Meanwhile, WL 12 (when 
the emission is detected, see variability in Fig.fTSTl and SVS 20 S 
show only a blueshifted component. In each case, the emission 
is relatively narrow, with FWHM between 10-50 km s . For 
most objects, the emission in this component is centered at the 
same spatial location on the detector as the continuum emission 
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Fig. 9. A comparison of the 12 CO P(6) and 13 CO R(6) line pro- 
files from GSS 30. The two profiles should be similar, includ- 
ing absorption from the same lower level. The absorption in the 
13 CO line is deeper and redder than that seen in the 12 CO line, 
despite that the 13 CO line cannot be more optically-thick than 
the 12 CO line. Instead, we infer the presence of redshifted ab- 
sorption from infalling gas at 10-20 km s . The broader emis- 
sion is produced far enough from the central source that it does 
not travel through the infalling gas and consequently does not 
suffer from absorption. 

and is generally not spatially extended. In one case, IRS 44 W, 
some of the 12 CO is offset by 0707 (~ 9 AU) W of the star and 
is spatially extended in the slit by Of.' 22 (~ 28 AU). Some very 
extended emission is also detected from GSS 30 and IRS 43 and 
is described in §3.3. In Appendix C, we discuss in detail the nar- 
row emission from GSS 30, IRS 44, and CrA IRS 2, which each 
have high S/N CO spectra. 

Figure [8] shows the rotational diagram for narrow emission 
from GSS 30, IRS 44 W, and CrA IRS 2. Relative to the broad 
emission component, the narrow emission component is pro- 
duced by gas that is cooler, with temperatures of 250-600 K, and 
more optically-thick, with column densities of log N(CO)~ 19. 
However, as indicated by the range in velocity centroids, the nar- 
row component may be produced by different processes for dif- 
ferent objects. For IRS 44 W and GSS 30, the temperature and 
column density from the rotational diagrams apply only to the 
blueshifted emission because 13 CO emission is not detected on 
the red side of the line profile. 

Figure compares the 12 CO P(6) fine with the 13 CO R(6) 
line for GSS 30. The two lines should have similar emission and 
absorption properties, modulated by the different optical depths. 
However, the 13 CO absorption is actually broader on the red side 
of the line than the 12 CO absorption. The most likely explanation 
for this discrepancy is that the CO emission suffers less from CO 
absorption than does the continuum. The stronger 12 CO emission 
then fills in the absorption more than the weaker 13 CO emission. 
The absorption, detected out to +10 km s , indicates the pres- 
ence of infalling gas in our line of sight to GSS 30. 

3.3. Spatially Extended CO Emission 

The M-band spectra of both GSS 30 and IRS 43 show CO emis- 
sion extended on arcsecond scales. We concentrate here on the 
extended CO emission from GSS 30 because it is much brighter 
than that from IRS 43. The slit position angle for the GSS 30 



observation roughly aligns wi th the ~ 45° position angle of the 
outflow (e.g. lAllen et"alll2002l) . 

Figure [10] shows that emission in low-/ I2 CO lines is de- 
tected out to 272 (265 AU at 120 pc) from the central source 
in both the NE and SW directions. Emission in 13 CO is also de- 
tected off-source in the SW direction and is barely detected to the 
NE. Co-adding C 18 0, CO v = 2-1, and CO v =1-0 high-/ lines 
each yields marginal detections of emission to the SW and non- 
detections to the NE. Most of the extracted emission is in narrow 
line profiles (FWHM ~ 7 km s _1 ) centered 2-3 km s _1 blueward 
of the CO absorption. A weak, broad (FWHM ~ 20 km s -1 ) base 
is also detected in the line profile. The CO absorption likely at- 
tenuates any emission on the red side of the 12 CO line profiles. 

Figure QT| compares the extent of emission in CO and the 
M-band continuum t o the spatial extent in archival L-band 
dDuchene et all l2007h and K-band images (3.6 and 2.2 /im, re- 
spectively) obtained with VL//NACO, deconvolved to the ap- 
proximate spatial resolution of our CRIRES observations. The 
K-band continuum and CO line emission is much more spatially 
extended than the L- and M-band continuum emission. The blue 
K-L color of the nebulosity suggests that the continuum emis- 
sion is seen in reflected light. 

Figure [8] shows the rotational diagram for CO line fluxes in 
the SW direction, extracted over a region from 0747 to 0793. 
Each continuum pixel in this region has ~ 2% the flux of the cen- 
tral pixel. A linear fit to the 13 CO emission yields a temperature 
of 250 ± 30 K. The low 12 CO/ 13 CO line flux ratio and the curved 
shape of the 12 CO rotational diagram indicates that the 12 CO 
emission is optically-thick, with logA^(CO) ~ 18.7. With these 
parameters, the C ls O is just below our detection limit. The line 
fluxes are reproduced with an emission area of (3.7 AU) 2 , which 
is two orders of magnitude smaller than the approximate extrac- 
tion region of 0743 x 073, or ~ (43 AU) 2 at 120 pc. Lowering the 
Doppler b parameter from 2.0 to 0.2 km s decreases the col- 
umn density to logA^(CO) ~ 17.8, thereby increasing the emis- 
sion area to (11 AU) 2 , which is still ten times smaller than the 
extraction region. This discrepancy between emitting area and 
extraction area is somewhat surprising because the CO emission 
appears to be smoothly distributed in the slit and the spatially- 
extended K-band continuum emission is smoothly distributed in 
the K-band image. 

3.4. Non-Detections of CO Emission 

Four objects, WL 6, IRS 44 E, Elias 32, and IRS 43 N, are un- 
detected in CO emission. For these four objects, the presence of 
strong CO emission relative to the continuum can be definitively 
ruled-out in our CRIRES spectra. However, the spectra are in- 
conclusive regarding the presence of weak CO emission. Each 
non-detection is discussed in detail in the following paragraphs. 

No CO emission is detected in our CRIRES spectrum of WL 
6. The telluric absorption falls on the red side of the line. In sev- 
eral other cases (WL 12, Elias 29, IRS 63), the CO v = 1 -0 emis- 
sion is detected only on the red side of the line because interstel- 
lar and wind absorption both occur in the same transitions. Some 
emission, with a strength that may be slightly above the contin- 
uum level, is detected between the telluric absorption and inter- 
stellar CO absorption. This weak emission could be either CO 
or an artifact from the telluric correction. If real, the weak CO 
emission would have blended with the interstellar/circumstellar 
absorption at low -resolution, leading to a non-detection in the 
ISAAC spectrum dPontoppidan et alfeoOSl) . The CRIRES spec- 
trum is consistent with weak CO emission. 
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Fig. 8. Plots of apparent column density versus excitation for CO emission for the broad and narrow emission components and 
spatially-extended emission. For the 12 CO v=l-0 lines, the synthetic fluxes are shown only for the P transitions. The dashed lines 
show the excitation expected from the N(CO) and T, based on models described in §A.4. 





Fig. 11. Top: VL77NACO K-band (left) and L-band (right, provided by Duchene) images of GSS 30, with the CRIRES slit in red. 
Nebulosity is detected to the NE and SW of the resolved central source. Bottom: The cross-dispersion profile of K-band, M-band, 
and narrow CO emission from GSS 30 along the CRIRES slit. The K-band emission is extracted from the NACO image, deconvolved 
to the 073 seeing during our CRIRES observation. The M-band continuum and the on-source CO emission are unresolved. That 
the M-band continuum emission is broader than the L-band continuum emission may be the result of broad wings on the CRIRES 
point-spread function, which we do not account for here. The extent of the narrow CO emission is somewhat similar to the K-band 
nebulosity. 



The non-detection of CO emission from IRS 44 E is compli- 
cated by the strong CO emission from IRS 44W. Weak emission 
is seen from the deep absorption to +50 km s~', with a peak-to- 
continuum that reaches ~ 0.05. However, this marginal detec- 
tion is not considered significant because it may be an artifact 
of the non-standard spectral extraction from this close binary 
(see §2.2). Regardless of whether this emission is attributable 
to IRS 44 E, most of the CO emission from the IRS 44 system 
is produced by the W component. Any CO emission from the 
E component is at least 5 times weaker than that from the W 
component. 

Bright emission was previousl y detected in high-7 CO l ines 
in an ISAAC spectrum of Elias 32 dPontoppidan et alll2003l) . so 
our non-detections of CO emission from Elias 32 on three dif- 



ferent nights 5 is surprising. Since CRIRES would spectrally re- 
solve CO emission from the narrow absorption line, the emission 
component should be detected at velocities that lack CO absorp- 
tion. The equivalent width required to explain the CO emission 
in the ISAAC spectrum is inconsistent with the non-detection of 
CO emission in the CRIRES spectrum of Elias 32. Either the 
equivalent width of the CO emission is variable or the ISAAC 
spectrum detected mostly extended CO emission. Some CO is 
detected in blueshifted absorption to Elias 32 (see §3.6). 



5 We observed Elias 32 three times, with each spectrum having low 
S/N. The highest S/N spectrum for Elias 32 was obtained on 3 May 
2008, when the telluric CO absorption was located at +23 km s~'. As 
with WL 6, the non-detection of CO emission on that date is not signif- 
icant. From the spectrum obtained on 10 August 2008, we find a limit 
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Fig. 10. Spectra around the 12 CO v = 1 - P(8) and 13 CO 
v = 1 — R(10) lines from GSS 30 seen on-source (black line), 
to the SW of the star (red line, multiplied by 5), and to the NE 
of the star (blue, multiplied by 5). The off-source spectra are ex- 
tracted from 0747 - 0793 from the central source. The 2D spec- 
tral images are obtained from coadding images over 12 CO and 
13 CO lines. Both CO and continuum emission are detected off- 
source in both directions along the slit. The emission is strongest 
to the SW. The off-source CO line emission is very narrow 
(FWHM~ 6.4 kms 1 ) and likely suffers from CO absorption 
in our line of sight. 

The object IRS 43 N is a very faint companion to IRS 43 
S. The non-detection is not significant, although the equivalent 
width must be at least 3 times smaller than the equivalent width 
of CO emission seen from IRS 43 S. 

3.5. H 2 S(9) Emission 

Emission in the H2 S(9) line at 4.6947 /im is detected from 9 
objects in our sample (see FigurefT2land Table|7]i. Most H2 lines 
have small or negligible blueshifts and a FWHM of ~ 20 - 30 
km s . Figure [13] shows the spatial distribution of the H2 emis- 
sion in the cross-dispersion direction, along with the spatial dis- 
tribution of continuum and CO emission. The H2 emission from 
five stars (IRS 63, Elias 29, SVS 20 S, L1551 IRS 5, and WL 6) 
is not extended along the slit direction and is consistent with a 

that the emission on the red side of the central absorption is less than 
7% above the continuum flux. Weak CO emission from several other 
embedded objects would not have been detected with this S/N. 



from Embedded Objects 



OS N 




-100 -50 50 100 
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Fig. 12. H2 line profiles from our sample, including non- 
detections from CrA IRS 2 and HL Tau. The H2 emission from 
SVS 20 S and WL 6 should be considered as tentative detections 
because the line is very broad. 



point source. The H2 emission is extended in the slit for the other 
4 objects with detected emission (GSS 30, WL 12, IRS 43, IRS 
44 W). The observed H2 emission from IRS 43 S is produced 
by molecular gas that extends > 1'.'5 (> 190 AU) to the S of the 
central object. The H2 emission from WL 12 extends 172 (~ 150 
AU) to the S of the object but is not seen to the N. 

Some of the H2 emission from IRS 44 is associated with IRS 
44 W, similar to the case for CO emission. Weak H2 emission is 
extended by up to ~ 171 (~ 140 AU) to the E of IRS 44 W. The 
equivalent width of H2 emission is ~ 50% smaller on our night 
of good seeing than on our two nights with poor seeing and than 
in the ISAAC spectrum (Table lE~3l . This result confirms that the 
S(9) emission from IRS 44 is spatially-extended because poor 
seeing increases the ratio of off-slit emission to on-slit emission. 

Some on-source non-detections are not significant. The CO 
3-2 R(10) line at 4.6958 /mi could mask any weak on-source H2 
emission from CrA IRS 2 and HH 100 IRS. The non-detection 
of H2 emis sion from E l ias 32 is limited by poor S/N. On the 
other hand. iBitner e t al. (2008) detected S(9) emission from HL 
Tau at a level that should have been detectable in our spectrum 6 . 
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Table 7. H 2 S(9) Line Detections 



Target 




FWHM 


EW 


Extended? 




km s -1 


km s _I 


km s _I 




Detected H 2 lines 


IRS 44 E+W 


-4.6(0.5) 


22(3) 


9.4(2.0) 


Y 


IRS 44 E+W c 


-4.9 (0.9) 


23(5) 


9.9(2.4) 


Y 


IRS 44 E+W 


-6.6 (0.6) 


25(3) 


6.7(1.3) 


Y 


Elias 29 


-9.5 (2.5) 


26 (8) 


0.9 (0.2) 


N' 


IRS 43 S+N 


-5.1 (1.6) 


18(4) 


1.7(0.3) 


Y 


IRS 63 f 


-10.2 (3.5) 


14(4) 


0.7(0.3) 


N 


WL I2 f 


-17(1) 


26(3) 


16(2) 


Y 


L1551 IRS 5 


-17(3) 


22(5) 


5.4(1.4) 


N" 


GSS 30* 


-7.9 (0.7) 


18(3.5) 


3.2(0.4) 


Y 


Broad H 2 lines: 


spurious detections?' 1 




WL 6 


1.5(2.0) 


89(10) 


6.9(1.7) 


N 


SVS 20 S 


19(2) 


70(14) 


5.3 (1.0) 


N 



Non-detections 



Elias 32 


low S/N 


CrA IRS 2 


blended with CO 


HH 100 IRS 


No detection, telluric correction? 


HLTau 


No detection 


RN0 91 


No detection 


TMC 1A 


No detection, bad telluric correction 


SVS 20 N 


low S/N 



"Relative to velocity of 1J CO and C 1!i O absorption. 
*2008-04-27 c 2008-04-30 ''2008-08-07 
Tow FWHM in continuum for CRIRES observation 
■flow S/N s C ls O line detected at -30 km s" 1 

''Uncertain detection because of poor telluric correction 



The H2 S(9) emission reported in iBitner et al.l d2008l) could be 
spatially-extended and not located within our slit, similar to the 
spati al distribution of H2 1-0 S(l) 2. 1 2 nm emission from HL 
Tau (iTakamiet all 120071: iBeck et al.L 120081) . However, in our 
2010 observation of HL Tau, the slit was aligned with the posi- 
tion angle of the outflow and the H2 1-0 S(l) 2.12 yu m emission, 
yet no extended H2 emission was detectable. The IBitner et"aT1 
(2008) spectrum of HL Tau includes a stronger unidentified fea- 
ture redward of the S(9) line, which may alternately point to a 
spurious detection in their spectrum. 

The broad, redshifted H2 emission from SVS 20 S is not con- 
firmed in the ISAAC spectrum and may be a spurious detection 
introduced by a poor telluric correction. A strong telluric CO2 
absorption line is located at +85 km s of the H2 line (+118 
km s" 1 relative to the source velocity at the time of the obser- 
vation). The H2 emission from WL 6 is also particularly broad, 
with a FWHM of 89 km s -1 . The ISAAC spectrum includes a 
similar detection for WL 6, so this broad emission may be real. 

The H2 S(8) line at 5.0529 //m is included in a wavelength 
setting that was used to observe only GSS 30. Any emission in 
this line is severely blended with strong emission in the CO v = 
1 - P(37) line. The H 2 S(8)/S(9) line flux ratio upper limit of 
~ 0.6 indicates that the emission is produced in gas hotter than 
~ 700 K. 

6 Using aFWHM=12 km s , as measur ed bv lBitner et aTJ ( 120081) . for 
a Gaussian line centered at the systemic velocity, we measure an equiv- 
alent width upper limit of ~ 0.1 km s -1 , roughly 2.5 times weaker than 
that detected in the TEXES spectrum. Some emission could be located 
at -10 to -30 km S" 1 in the on-source spectrum but not be detectable 
because of a poor telluric correction. Our March 2010 observation of 
HL Tau has a poor telluric correction at the H 2 line, presumably be- 
cause of variable telluric absorption and is not usable for the analysis 
of on-source emission. No off-source H 2 emission was detected in that 
observation. 




Fig. 13. Spatial profiles of M-band continuum emission, coadded 
12 CO v = 1 — line emission, and H2 S(9) emission from four 
sources. The H2 emission from these four sources is much more 
extended than the CO and continuum emission, an d is reminis- 
cent o f H2 v = 1-0 S(l) emission from outflows dBeck et al.L 
120081) . The H 2 emission from Elias 29, L1551 IRS 5, IRS 63, 
and WL 6 is not extended and not shown here. The CO emission 
from GSS 30 shown here is the blueshifted component between 
-20 and -10 km s _1 . The CO emission on the line wings is not 
spatially extended beyond the continuum emission. 



3.6. CO Wind Absorption 

Within our sample, six objects (HH 100 IRS, IRS 63, WL 12, 
Elias 29, Elias 32, and TMC 1A) show wind absorption in CO 
v = 1 - lines (see Figs, g] and [H). The wind to HH 100 IRS, 
TMC 1A, and Elias 29 is seen to velocities of ~ 100 km s . 
The wind absorption to Elias 32, IRS 63, and WL 12 is slower, 
with speeds up to ~ 50 km s _1 . These differences, and the lack 
of CO wind absorption in the majority of our sample, may be 
attributable to the different viewing inclinations, compositions, 
and wind velocities between sources. 

The wind absorption is typically weakest at J" = and peaks 
in strength at J"= 10 to 20. For HH 100 IRS, the absorption is 
strongest in two distinct components at -65 to -85 km s _1 . The 
absorption seen to TMC 1A is strongest at low velocity and be- 
comes undetectable at < -65 km s _1 . To calculate temperatures 
and column densities, we assume that the wind absorption is op- 
tically thin. An alternate explanation, which applies to the slow 
wind seen from TMC 1 A described below, is that the wind only 
attenuates a fraction of the M-band continuum emission. In this 
case, the absorption in low J and high J lines would still need 
to be optically thin because the absorption depth is largest for 
lines with 10 < J < 20. For TMC 1A, an optically-thick, low- 
velocity (< 20km s _1 ) component is seen in low-/ lines and in 
I3 CO lines. Since the spectrum does not go to at low velocities 
in the 12 CO lines, the low-velocity absorption must obscure only 
about half of the M-band continuum emission produced by the 
inner disk. 

The temperature of the wind for HH 100 IRS and TMC 1 A 
is 1260 + 100 K and 1260 + 40 K, respectively (Fig. [B). Fig. [BJ 
shows the column density of CO in 1 km s bins in the winds 
seen to HH 100 IRS and TMC 1 A. The log of the total CO col- 
umn density to TMC 1A, summed from 20-65 km s , is 19.6, 
and for HH 100 IRS, summed from 60-90 km s 1 , is 18.4. 

FigurefTBIshows CO v = 1-0 emission from observations of 
HH 100 IRS obtained 11 months apart, and of WL 12 obtained 
3 years apart. For HH 100 IRS the blueshifted CO absorption 
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Fig. 15. Rotational diagram for wind absorption from HH 100 IRS (left) and TMC 1 A (middle) covering the -55 to -80 km s _1 
region that includes the double-peaked wind absorption. The temperature of the absorbing gas for both stars is 1260 + 100 K. The 
right panel shows the column density of CO in 1 km s _1 bins. 



4. DISCUSSION 
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Fig. 14. Wind absorption towards four sources in 12 CO 1-1 lines, 
with high- to low-/ from top to bottom. 



changes between the two observations, with stronger absorption 

at 70 km s in August 2008. The shape and strength of the 

emission line profile remained the same, to within observational 
uncertainties. For WL 12, some of the absorbing gas appears to 
have gone into emission, which could happen with a change in 
viewing angle of the wind or if the opacity of the wind decreased, 
allowing some emission to escape. 



In §3, we describe four different components of CO ro vibra- 
tional lines from embedded young stars: (1) broad, vibrationally- 
excited emission (§3.1), (2) very narrow, spatially-extended 
emission (§3.3), (3) absorption from winds and envelopes (§3.6), 
and (4) narrow, optic ally- thick emission (§3.2). In the following 
subsections, we discuss the physical interpretations of the differ- 
ent components of warm CO gas in the YSO 

4.1. Disks in Broad CO Emission 

The broad CO emission seen from the embedded objects in our 
sample is similar to the warm CO emission seen from CTTSs. 
Of the 12 CTTSs with detected CO v = 1 - emission in the 
Naiit a et alj(l2003l) sample, nine are also detected in CO v = 2— 1 
lines. For these CTTSs, the FWHM range from 56-134 km s _1 
and the line fluxes range from (0.5 - 14) x 10~ 14 erg crrT 2 s _1 
(0.3 - 10 x 10~ 5 L Q at 140 pc). Similarly, in our sample most ob- 
jects with broad CO v = 1 — line emission also show broad CO 
v = 2 — 1 line emission, with FWHM that range from 55-140 
km s _1 and line fluxes that range from ~ (2 — 10) x 10~ 14 erg 
crrT 2 s _1 (2 - 10 x 10~ 5 L for a 120 pc distance and ext inction 
A M = 1 mag.). In the CTTS sample of lNaiita et all (120031) . I3 CO 
emission was only detected to one source. In our spectra, the 
broad component is not detected in 13 CO lines. The measured 
temperatures are roughl y consistent with the CO temperatures 
in the CTTS sample in lSalvket all (120091) . The centroids of the 
spectral line profile are consistent with the source velocities, and 
the emission is centered on the source and not spatially extended 
beyond ~ 10 AU. The most likely origin for this broad compo- 
nent of CO emission from embedded sources is the disk. 

That vibrational levels v = 2,3,4 are populated despite rota- 
tional temperatures of ~ 1000 K supports the disk interpretation 
for this emission, since large populations i n the v > 2 leve l of CO 
typically requires FUV photoexcitation dKrotkov et al.L 1 1980; 
Ivan Dishoeck & BlackL 119881: iBrittain et all 120071) . Electronic 
(A-X) transitions of CO in the FUV are much stronger than the 
rovibrational transitions and become optically-thick at smaller 
optical depths. As a consequence, gas that produces UV-excited 
emission from v' > 2 will not be detectable in 13 CO. The lack 
of broad 13 CO emission is therefore consistent with the presence 
of UV-excited gas and places a limit of logA^(CO)< 18 in this 
inner disk region. 

Emission from the v' = 2 level has previously been detected 
from emb edded young st ars in v = 2-0 overtone transitions from 
WL 16 (ICarr et al.L [l993). in overtone and fundamental emission 
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Velocity (km s ) Velocity (km s ) 

Fig. 16. Left: Comparing the CO fundamental line profile from HH 100 as seen in observations separated by 1 1 months. The profiles 
shown here are the co-added lines of P(4), P(6), P(7), and P(8). No significant difference is seen in the emission line profile between 
the two epocsh. In August 2008 the wind absorption at -75 km s _1 was deeper than that seen in September 2007. Right: Co-added 
CO emission from two observations of WL 12. The observation of September 2007 covered low-J lines while the observation of 
March 2010 covered high-J lines, so the comparison of line profiles is not direct. Some of the wind absorption in 2007 appears to 
go into emission in 2010. 



from CTTS s dNaiita et all l2003t iBast et al.L I20TTI) and Herbi g 
AeBe stars (iBrittain et al 1 l2007l l2009h Ivan der Plas et aU l2009). 
That emission from v' > 1 traces the broad component is con- 
sistent with the far-UV field being strongest near the star. In 
contrast, far-UV-excited CO rovibrational emission from Herbig 
AeBe stars is usually produced at larger disk radii. UV radiation 
fields from embedded young stars are not directly detectable but 
should be sufficiently strong to excite warm CO in the inner disk. 



Figure [17] compares the rotational diagram of number of 
CO molecules for HH 100 IRS, DR Tau, and TW Hya . The 
fluxes for TW Hya were measured bv lSalvket all (120071) . The 
fluxes from DR Tau were measur ed from the CRIRES spec- 
trum presented in [Bast et al] d2011l) . TW Hya is an ~ 8 Myr star 
(Mam aiel 120051) with a disk that has an inner dust hole of ~ 4 
AU dCalvet et al.L |2002|) and an accretio n rate of ~ 1.5 x 10~ 9 
M Q vr" 1 dHerczeg & Hillenbrand! l2008h . DR Tau is a younger 
star w ith an accretion rate of 3 x 10 -7 M Q yr" 1 dGullbring et all 
2000), comparable to the 10~ 6 M yr 1 estimated for HH 100 
IRS bv lNisini et all d2005l) . The strength of CO v = 1-0 emis- 
sion from DR Tau is remarkably similar to that from HH 100 
IRS. The inner disk structure of HH 100 IRS is probably some- 
what similar to that of DR Tau and other CTTSs with high ac- 
cretion rates. The CO emission from both HH 100 IRS and DR 
Tau is much stronger than that from TW Hya, perhaps because 
TW Hya has a smaller accretion rate and therefore less viscous 
heating in the inner disk. 



Double-peaked profiles are absent within our sample and are 
also somewhat rare for CTTSs. Of the many previously pub- 
lished samples of CO emission from CTTSs, only AA Tau, SR 
21, GQ Lup, RNO 90, Elias 23 (this work, see Fig. 5), and DF 



Tau are resolved as double-peaked profiles 7 ( Naiita et al 
Salvk et all 120071; IPontoppid an et all l2008HSalvk et al 
Naiita et all 120091: iHugelmev er et all l2009t IBast et al. 



2003 



2009; 



2011 



Brown et all lin prep.1) . Most CO lines from embedded objects 
are consistent with Gaussian profiles. Some Gaussian profiles 
may be consistent with a Keplerian disk viewed face-on, as is 



the case for CO lines f rom several CTTSs (Pontop pidan et all 
120081 iGoto eTaTll201 ll) . 

To investiga te the absence of double-peaked line profiles, 
IBast et all d201 lb selected a subsample of 8 CTTSs with bright 
CO lines and high accretion rates. In that sample, the combi- 
nation of the narrow central peak and the lack of any spatially- 
extended emission, with upper limits of a few AU, indicates the 
presence of some gas near the star that is not in Keplerian rota- 
tion. The subset of embedded objects that have low bolometric 
luminosities is likely similar to the sample of Bast et al.. In both, 
the broad line widths indicate an origin near the where the in- 
ner disk is li kely trunca t ed. Th e fluxes suggest an emitting area 
of ~ 1 AU 2 . IBast et all j20TTh were unable to explain the line 
fluxes, profiles, and lack of any spatially-extended emission with 
an origin in the surface layers of a warm disk in Keplerian rota- 
tion and with a surface temperature described by a power law. 
One possi ble explanation attributes much of the line flux to a 
disk wind dPontoppidan et all 1201 lb . 

Regardless of explanation, the large velocities in the line pro- 
file are likely explained by Keplerian broadening. For the disk 
explanation, the Keplerian broadening is in situ, while for the 
disk wind explanation the Keplerian broadening would apply 
to the inner radius of the launch region. Within the context of 
the disk/disk wind interpretation, the inner radii /?,„(sin i) 2 for 
the CO emission are calculated from the velocity of the best-fit 
Gaussian profile at 20% the peak flux (Table [8]) and a central 
mass 1 M e . Although a direct measurement of this location in 
the line profile is typically not possible in our sample because 
of insufficient S/N, all broad line profiles are consistent with a 
Gaussian profile. About 93% of the warm CO emission is pro- 
duced beyond this inner radius. Alternately, we solve for inclina- 
tion by setting Rj„ = 0.05 AU for every source. The listed inner 
radii and inclinations suffer from significant uncertainties. 



7 iNaiita et"al] (2003) suggest that the CO emission from RW Aur is 
double-peaked, but this inference was based on modeling a small spec- 
tral region that did not show a clear double-peaked line profile. We 
find no evidence for a double-peaked CO emission profile in a CRIRES 
spectrum of RW Aur A (Brown et al., in preparation). 
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Fig. 17. A comparison of 12 CO emission from the broad compo- 
nent of HH 100 IRS (black circles), DR Tau (blue diamonds), 
and TW Hya (red squares). The emission in 12 CO v = 1 — 
is very similar for HH 100 IRS and the CTTS DR Tau, with 
much weaker emission seen from the older CTTS TW Hya. CO 
V = 2 — 1 emission is stronger from HH 100 IRS than DR Tau. 



4.2. Molecules in Winds H 

Mole cular emission is o f ten detected in winds from young star s 
(e.g. IShang et all l2007t iBeck et all 120081: iDavis et all l2010h . 
As described in IPanoglou et all dsubmittedl) . winds can have a 
molecular component for one of the following reasons: (i) en- 
trainment of molecular gas in the envelope or cloud by the wind, 
(ii) molecular formation within the wind itself, or (iii) launch- 
ing of molecular gas in a disk wind. Typically, observations of 
molecules in winds lack sufficient spatial resolution to discrim- 
inate between these possibilities. Our data includes diagnostics 
of the wind near the launch region and also of the slow-moving 
molecular gas extended from the star. 



4.2.1 . Molecules launched in the wind 

The fast blueshifted CO absorption detected to HH 100 IRS, 
Elias 29, and TMC 1A suggests that the MHD wind itself is 
at least partially molecular when it is launched. The slower CO 
absorption detected to Elias 32, WL 12, and IRS 63 could be 
explained by either an MHD disk wind or a slower thermal disk 
wind. When detected, wind absorption usually occurs close to 
the star because the density within the wind decreases with ra- 
dius squared as the wind expands. Molecular gas from the enve- 
lope that gets entrained in the outflow is highly unlikely to pro- 
duce absorption close to the star and at such high velocities. CO 
wind absorptio n in the M-band has also been detected from high- 
mass YSOs dMitchell et all 1 19901: iThi et all l2010h . The pres- 
ence of CO in the wind suggests that the wind from these sources 
favors launching fr om the disk ra t her than the stellar chromo - 
sphere or corona dMatt & Pudritzl 120051: iFerreira et all [2006). 
Although any significant ma ss loss from the star i tself would 
have to be cooler than 10 4 K ( Mat t & Pudritzl l2007l> . a chromo- 
spheric wind would likely be mostly neutral. A s an alternative 
to a w ind that is molecular where it is launched, Glassgold et all 
(1199 11) suggest that CO, but not H2, may form within the wind 



itself. However, in their models a high CO abundance requires 
a mass loss rate > 10~ 5 M yr , which is lar ger than the mass 
loss r ates inferred for typical Stage 1 sources dBontemps et all 
1 19961) . 

Models of MHD disk winds by Panoglou et al. ( submitted) 
demonstrate that molecules can survive within the wind when 
the mass loss/accretion rate are sufficiently large (M acc ~ 10~ 6 
M e yr -1 ), in material with an ionization fraction that is high 
enough to couple the molecular gas to the magnetic field. The 
measured temperature of 1260 ± 100 K is roughly consistent 
with the predicted wind temperatures for objects with accre- 
tion rates of 10~ 7 - 10~ 6 M yr -1 . Wind absorption in CO has 
not been detected previously to CTTSs. Several CTTSs show 
on-sour ce FUV H2 emission with velocities up to 50 - 100 
km s _1 dHerczeg et all 12006). The frequency of CO absorption 
in winds from embedded objects and deficiency of CO absorp- 
tion in winds from CTTSs is consistent with the survival of CO 
requiring large outflow rates to shield the outflow from irradia- 
tion by the central star. 

That CO wind absorption is detected toward only some of 
the objects in our sample is consistent with the MHD wind 
being non-spherical and with a relatively wide opening angle 
near the star. In a sur vey of the He I /U0830 line from CTTSs, 
Edwards et al ] (I2006h found that winds from stars seen pole-on 
tend to be faster and more optically-thick than winds from stars 
seen closer to edge-on. Within our sample, CO wind absorption 
is not detected to the two stars with the broadest CO emission 
lines, indicating that perhaps the same relation holds. On the 
other hand, the narrow CO line emission from CrA IRS 2 in- 
dicates a low disk inclination, but no CO absorption is detected. 
The lack of CO wind a bsorpt i on from HL Tau, with an incli- 
nation of 65-70° dClose eta Il997t iLucas et all 12004ft 9 . also 
contr asts sharply w i th the deep, fast He I absorption seen to the 
star dEdwards et all 120061) . For both CrA IRS 2 and HL Tau, the 
likely interpretation is that the molecular fraction is low in the 
wind. For HL Tau, the discrepancy between the deep He I wind 
absorption and the lack of any detectable CO wind absorption 
could also be explained if 4.8 /im continuum is seen through a 
very different line-of-sight than the 1 . 1 fim continuum emission. 



4.2.2. Spatially-extended molecular emission 

Powerful jets and winds from young star s can carve out cavities 
within the circumstellar envelope (e . g ., Whitney & Hartmann. 



119931: IWoodet all 1200 UlYbarraet all r2006). The interaction re- 
gion between the cavity and the jet/wind entrains some cold gas, 
which is likely seen in bipolar molecular outflows with veloci- 
ties of a few km s _1 . At the surface of the cavity wall, gas can 
be heated by shocks and UV photons from t he central star (e.g. 
ISpaans et all 1 19951: Ivan Kempen et alll2009h . 

Within our sample, H2 0-0 S(9) emission is detected in nine 
objects. For four of these detections the H2 emission is spatially- 
extended, always asymetricalry about the star. The spectral line 
profiles are typically centered within 10 km s _1 of the systemic 
velocity. The extended component is likely produced in the en- 
velope or nearby cloud material that is shocked by powerful out- 



8 For simplicity, the term winds used here is intended to incorporate 
all disk winds, MHD winds, and jets from a YSO. 



9 dFurlan et all 12008) suggest an i ~ 15° for the HL Tau disk from 
modelling the SED, but we consider the near-IR polarimetry a more 
reliable method to measure inclination because inclination is degenerate 
with many other parameters in broadband SED fitting and because the 
combination of high extinction and lack of ice absorption suggests disk 
attenuation rather than envelope extinction. 
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flows. iGreene et all (120101) detected H 2 1-0 S(0) and 1-0 S(l) 
emission from all 17 embedded young objects in their sample, of 
which 10 have emission that is that is spatially extended on 2-3" 
scales. Compared with the H2 0-0 S(9) line reported here, the H2 
1-0 S(0) and 1-0 S(l) emission lines have similar FWHM (20- 
30 km s ) and central velocities that differ by ~ 5 - 10 km s _1 . 

Of the objects with extended H2 emission, the slit PA was 
well-aligned (better than 20°) with the outflow axis only for 
GSS 30. The spatial distribution is likely analogous to the 
distribution of warm H2 emission seen towards other embed- 
ded YSOs and e nvironmentally-young CTTSs that drive power- 
ful outflows (e.g . McCaughrean & Ma c LowLll997LlDa vis et al., 
200UIWalter et aUl2003l:lBeck et al l 120081; iNeufeld et aUl2008 ; 
Lahui s et al Greene et al .. 2010). The spatially unre- 



Table 8. Inner Radius of disk CO Emission 



solved H2 emission could be produced by a disk, as is the 
case for w e ak H 2 rovibrational emission from a few CTTSs 
dBarv et all 12008?). However, because the H2 emission of- 
ten includes a significant contribution from surrounding enve- 
lope/cloud material, a careful analysis would be required to use 
the line as a disk probe. In light of our results, the origin of H2 
emission from GS S 30, and perhaps Elias 29 and HL Tau, by 
Bitn er et al.l (|2008) should be considered the wind/envelope in- 
teraction region rather than the disk. 

Spatially-extended CO emission is also detected from two 
objects, GSS 30 and IRS 43, that show extended H2 emission. 
The extended CO line emission traces cooler gas and is spec- 
trally more narrow than the H2 line emission. That extended CO 
emission is not detected more frequently is likely explained by 
densities that are much lower than the critical density of 5 x 10 12 
cm -3 required to populat e the v _ 1 levels n eeded to produce CO 
rovibrational emission dNaiitaet all 1 1996b . although instead the 
CO/H2 abundance ratio could be negligible. On the other hand, 
the presence of extended CO emission from two objects may in- 
dicate high densities in the associated outflows. 

4.3. Disks and/or Outflows in Narrow CO Emission? 

Narrow CO emission is detected from 9 of 18 embedded ob- 
jects within our sample. The optical depth of this component 
is larger than that in the broad component, as evidenced by the 
high 13 CO/ 12 CO line ratios. Narrow line profiles have a wide 
range of properties, indicating that this classification is overly 
broad. For GSS 30, IRS 44 W, WL 12, and SVS 20 S, some or 
all of the narrow CO emission is blueshifted by ~ 10 km s _1 . 
For other objects the emission centroid is consistent with the 
systemic velocity. The lack of spatially-extended emission indi- 
cates that this emission is produced relatively close to the central 
star, although some narrow emission from IRS 44 W is slightly 
extended, and 12 CO emission from GSS 30 is likely emitted be- 
yond infalling CO absorption. The different velocities of the nar- 
row component indicates the emission is produced by different 
physical processes, depending on the star. 

This narrow componen t does not have a counterp art for the 
typical C TTSs observed bv lNaiita etail d2003l) and lBrown et alj 
din prep.1) . The mechanism that produces an optically-thick col- 
umn of ~ 400 K gas for these evolutionary young objects does 
not heat a similar amount of gas in more evolved CTTSs. On 
the other hand, 13 CO emission is commonly detected within the 
iBast et ail d201 ll) sample of high accretion rate CTTS. In par- 
ticular, the CO emission from S CrA S is remarkably similar to 
that seen from HH 100 IRS, IRS 63, and TMC 1 A. All four ob- 
jects show broad emission in CO v = 2 — 1 transitions, narrow 
emission in 13 CO transitions, and optically-thick 12 CO v = 1 — 
transitions that include both components. The CTTSs with large 
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accretion rates are likely not much more evolved than many of 
the embedded sources in our sample, including HH 100 IRS, 
TMC 1A, and HL Tau. Four Stage I YSOs with high luminosi- 
ties show some blueshifted CO emission, which indicates that 
the mechanism that produces the warm CO continues to change 
with increasing bolometric luminosity. 

The component could be explained by a thick layer of gas 
that is either in a thermal disk wind and/or the disk surface. 
Evaporative winds typic ally produce emission lines w i th cen - 
troids at - 10 km s ' dAlexanden 1 120081: lOwen et all l2010h . 
Different viewing anlges could explain the range of measured 
velocity centroids. The lack of significant emission in lines from 
v' > 2 indicates that the CO emission region is shielded from the 
UV emission from the central star, which should be bright. 

4.4. What Determines the Strength of CO Emission from the 
Disk? 

Most of the objects with large bolometric luminosities emit in 
the narrow component but not in the broad component (see 
Tableland Fig. [T8l . The sources with small bolometric lumi- 
nosities typically emit in the broad component. A few objects 
with intermediate bolometric luminosities emit in both compo- 
nents. Within these categories, the correlation between bolomet- 
ric luminosity and the strength of the broad component (v = 2 - 1 
line flux) is poor, while the bolometric luminosity may be corre- 
lated with the strength of the narrow component ( 13 CO line flux). 
These comparisons are not definitive because the sample size is 
small and the bolometric luminosities have large uncertainties. 
Other YSO properties must also determine the strength of the 
emission, as IRS 44 E lacks strong CO emission despite a high 
bolometric luminosity. However, bolometric luminosity clearly 
affects what regions emit in CO fundamental transitions. 

For CTTSs and low-luminosity embedded objects, the broad, 
vibrationally-excited emission indicates that molecules survive 
in the inner disk near the disk truncation radius. On the other 
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Fig. 18. The luminosity in the I3 CO v = 1-0 R(13) line and 
the CO v = 2 — 1 R(7) line versus bolometric luminosity. The 
13 CO and v = 2 — 1 lines trace emission in the narrow and broad 
component, respectively. Upper limits were calculated assuming 
a FWHM of 30 km s _1 in the narrow component and 100 km s _I 
in the broad component. 



hand, the non-detection of this broad emission from some 
embedded objects indicates an absence of warm CO gas in 
Keplerian rotation near the star. Of the high-luminosity objects 
in our sample that lack broad CO emission lines, disks have 
been detected around GSS 30 and L1551 I RS 5 in the sub-mm 
(lLoonev et allll997tlJ0rgensen et all 12009) and around GSS 30 
Elias 29, and IRS 44 E in near-IR polarimetry dBeckford et all 

2008) . The measured di sk inclinations of Eli as 29 (~ 53°) and 
GSS 30 (~ 66°) from iBeckford et all d2008l) rule out that the 
narrow emission simply results from disks seen face-on. 

For YSOs with high bolometric luminosity, molecules may 
be photodissociated out to larger radii than is typical for CTTSs. 
As a consequence, perhaps no warm, broad CO emission is 
detected because the inner disk is mostly atomic or ionized. 
Similarly, Herbig AeBe stars frequently emit in CO fundamental 
lines at large radii (e.g. Brit tain et all 120091: Ivan der Plas et all 

2009) but in trac ers of warmer ga s, such as [O I] /16300, from 
their inner disks dAcke et alll2005l) . As with some Herbig AeBe 
stars, the line widths of the narrow component from Stage I 
YSOs are consistent with emission from larger radii in the disk. 
A search for broad emission in tracers of neutral and ionized gas 
could test whether the disks around these stars simply have a dif- 
ferent temperature structure than that of CTTSs or whether the 
velocity structure of the young, forming disk differs from that of 
mature disks. For stars with blueshifted CO emission, the region 
where the CO emission is produced may be active or turbulent 
enough so that the disk surface evaporates. The deficiency of 
strong, UV-excited lines from v' > 1 for high bolometric YSOs 
indicates that the warm CO gas is not irradiated by bright UV 
radiation. The bright UV emission expected from central source 
would need to be attenuated, perhaps by dust in the wind, before 
reaching the CO gas. 

While the inner disks around low-luminosity embedded ob- 
jects are likely similar to those around CTTSs, the inner part 
of any disks around embedded objects with moderate luminosi- 
ties (higher mass and/or younger) are probably quite different. 
Any CO gas in the inner disks of these objects should have been 
detected. If the inner disk has not yet settled into Keplerian ro- 
tation, then CO gas in the inner disk, if present, could produce 



the narrower, optically-thick emission. Alternately, dust that is 
optically-thick at 4.8 fim could shield any warm CO emission, 
particularly if the inner disk does not have a temperature inver- 
sion because of significant viscous heating. 

5. CONCLUSIONS 

We present an overview of fundamental CO emission in high- 
resolution VLT/CRIRES spectra of 18 embedded young stars, 
obtained to explore the warm inner disks as they are still form- 
ing. CO emission was detected from 14 of the 18 stars. The four 
non-detections are only significant in ruling out large equivalent 
widths in the lines. We find the following results: 

1 . On-source CO emission is produced in narrow and broad 
components, both of which are single-peaked. Often only one of 
these components is detected from a given source. The narrow 
emission component is typically detected from embedded ob- 
jects with high bolometric luminosity, while the broad emission 
component is typically detected from objects with low bolomet- 
ric luminosity. 

2. The broad emission (FWHM range of 50-160 km s _1 ) is 
detected in lines v = 1— 0, v = 2 - I and, when the S/N is suf- 
ficient, lines with v' > 2. This emission traces warm (~ 1000 
K) gas covering ~ 1 AU 2 and is centered at the systemic veloc- 
ity. This emission is reminiscent of the warm, single-peaked CO 
emission detected from more mature disks around CTTSs with 
high accretion rates, and likely traces gas in the inner regions of 
the disk or in a slow, thermal disk wind. The vibrational excita- 
tion is higher than the rotational excitation, which suggests that 
UV photoexcitation of the disk surface may be responsible for 
emission from levels with V > I. 

3. Narrow emission (FWHM range of 10-50 km s" 1 ) is de- 
tected in 12 CO, 13 CO and, when S/N is sufficient, C ls O and CO 
v = 2 — 1 lines, indicating large optical depths. This emission 
tends to be cool (300^1-00 K) and is usually not spatially re- 
solved, except in the case of IRS 44. Some of these lines are 
blueshifted by 5-15 km s , indicating an origin in a slow out- 
flow. Other lines are at the systemic velocity, either because the 
gas is not escaping the disk or because of viewing angle of the 
wind. The range of line profiles indicate that this classification 
includes multiple physical processes. The faint emission in CO 
v — 2 - 1 lines indicates that this CO emission is shielded from 
bright UV emission produced by the central star. 

4. Broad CO emission is not detected from objects with bolo- 
metric luminosities > 15 L Q . Either warm CO gas is not present 
in the inner disk of these objects due to photodissociation, any 
warm CO is shielded by dust that is optically-thick at 4.8 fim, or 
the disk has not yet settled into Keplerian rotation. If the molec- 
ular fraction is low near the star, the gaseous disk should be de- 
tectable in atomic or ionized gas lines. 

5. Very narrow, spatially extended CO emission from GSS 
30 and IRS 43 likely traces cold molecular gas, perhaps in dense 
regions where the winds interact with the envelope and nearby 
molecular material. The extended CO emission is especially 
bright from GSS 30, indicating some unique physics or mor- 
phology for that source which is not typical of the rest of our 
sample. 

6. Winds are detected in CO absorption towards 6 ob- 
jects, with velocities of 10 to 100 km s and temperatures of 
1260 + 100 K. Because wind absorption usually occurs close 
to the star, these detections indicate that the wind is partially 
molecular when it is launched. A molecular wind is consistent 
with a wind launched from the disk but not with an accretion- 
driven wind launched from the stellar chromosphere or corona. 
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The fraction of objects with CO detected in the wind (6 of 18) 
suggests that the wind may have a wide opening angle near the 
star. 

7. H2 S(9) emission is commonly detected in our sample. The 
spatial extents of ~ 1 - 2" seen in about half of these detections 
indicate that the H2 emission is produced in winds interacting 
with surrounding material rather than the disk. 

We have analyzed M-band CO emission as a probe of gas 
in the inner few AU of disks around young stars that are still in 
the embedded phase of pre-main sequence evolution. The data 
quality in these observations, obtained with a sensitive, high- 
resolution instrument on an 8m telescope on nights with good 
seeing, will be difficult to improve upon with existing instru- 
mentation, although larger samples and spectral imaging with 
current capabilities could clarify the interpretation of narrow CO 
emission and spatially-extended CO emission. Spectral imaging 
with future instruments, such as NIRSpec on JWST or METIS 
on the Extremely Large Telescope, holds significant potential 
for understanding the morphology and production mechanism 
of CO emission around GSS 30 and other sources, while spec- 
troastrometry of embedded objects obtained with a laser-guided 
AO system could identify th e location of the different CO emis- 
sion components, following Pontoppi dan et al.1 (120081) . Despite 
its revolutionary capabilities in observations of young disks, 
ALMA will be unable to trace warm gas within a few AU of the 
central star, where terrestrial planets are thought to form. 
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Appendix A: Resolved Binaries 

In our long-slit spectra, spatial information is obtained in the 
cross-dispersion direction. Most observations are consistent with 
point-source emission. The objects IRS 44, IRS 43, and SVS 20 
are known binaries that were resolved in our data and are de- 
scribed below. Interferometry of LI 551 IRS 5 at 7-mm indicates 
two di stinct sources of dust continuum emission separated by 
~ 073 (lLoonev et allll997LlRodriguez et al.Lll998l) . although the 
two components were not resolved in hig h-resolution H, K, or L- 
band imaging by (Duch ene et al.L f2007). Regardless of whether 
L1551 IRS 5 is a binary, the continuum has a FWHM of 078 in 
our data (either poor seeing or spatially-extended M-band con- 
tinuum emission), so resolving the two components is not possi- 
ble. 



A.1. IRS 44 

IRS 44 w as definitively identifi ed as a close binary in L-band AO 
imaging dDuchene et al.Ll2007h . which revealed two components 
with a flux ratio close to 1.0 and a separation of 0730 at a PA 
of 87°. High-resolution K-band images also suggested that IRS 
44 is a binary with a similar separation and a flux ratio of ~ 5 
(Ratzka et al. 2005; see also Terebey et al. 2001, Allen et al. 
2002). However, the K-band imaging is somewhat less reliable 
than the L-band imaging because the large binary flux ratio and 
bright nebulosity leads to poor S/N and uncertain detections of 
the secondary. 

Of our five CRIRES observations of IRS 44, the binary was 
resolved only on the nights of 6-7 Aug. 2008, when the see- 
ing was superb (~ 0732). The two components are marginally 
resolved with a separation of about 073 on the detector (see 
Fig- US. This separation does not include the ~ 13° offset in po- 
sition ang le between the slit an d the binary. Assuming the binary 
PA from iDuchene et al.1 (120071) . then IRS 44 E is ~ 0.69 + 0.12 
mag brighter than IRS 44 W in the M-band continuum, includ- 
ing a 0.08 mag uncertainty attributed to uncertainty in t he slit 
position relative to the two components. Duche ne" et al.l (120071) 
suggested based on the K- and L-band observations that IRS 44 
W was more heavily reddened than IRS 44 E and is the primary 
in the system. This scenario predicts that the IRS 44 W would 
be the brighter component in the M-band, which is the opposite 
of what is seen. 

In the CO 4.67 pm ice band, absorption reduces the con- 
tinuum flux from IRS 44 E by ~ 18% and from IRS 44 W by 
~ 22%. Some ice absorption must be local to IRS 44 W. 

A.2. IRS 43 

IRS 43 was resolved as a 0757 bi nary with an L-band fl ux ratio 
of 3.1 mag and a PA of 336° bv lDuchene et al.1 (120071) . In our 
M-band long-slit spectra, the continuum emission from IRS 43 
is concentrated on the primary. Some extended emission is de- 
tected to the N, in the direction of the faint secondary. If this 
emission is attributed to the secondary component, then it is 
~ 1.8 + 0.5 mag fainter than the primary in the M-band con- 
tinuum. 

The 4.67 pm CO ice band absorbs 80% of the photons emit- 
ted from both the continuum and from the secondary (or ex- 
tended nebulosity) to the N. We infer that the ice absorption 
likely occurs in a cloud that envelops both objects. 

A.3. SVS 20 

SVS 20 is a 1751 binary with a PA of 9.9° (lEiroa et al.L \\9S% 
lHaisch et all 120021 l2006h . The separation on our slit is ~ 1756, 
with an M-band continuum magnitude difference of 1 .25, similar 
to the magnitude difference seen at other optical and IR wave- 
lengths. The luminosity ratio of ~ 500 listed in Table 2 is likely 
much too high. 

The properties of the stellar components of SVS 20 are un- 
certain. SVS 20 S is considered the primary component be- 
cause it is brighter in the near-IR and mid-IR. Weak K-band 
absorpt ion features indicate th at the primary is likely an early- 
G star dDoppmann et all 120051) . although as with IRS 63, this 
spectral type is tentative because the photospheric velocity is 
discrepant by 10 km s - 1 from the v/ y of th e HCO + emission 
( Greg ersen et all l200dh . lOliveira et al.1 (120091) found that one of 
the components is an M4 star with Ay — 3 mag., based on TiO 
features in an optical spectrum that included unresolved light 
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from both objects. The low extinction suggests that whichever 
component is the M-dwarf is not embedded in the envelope and 
perhaps dominates the optical flux as a result. The 13 CO ab- 
sorption lines are also more optically thick to SVS 20 S than 
to SVS 20 N. On the other hand, the CO ice absorption has a 
larger opacity, r ~ 1, to SVS20 S than to SVS 20 N (see also 
Pontoppid an et alll2003l) . The brightness difference between the 
two components, ~ 1.3 mag in the near-IR and ~ 1.0 mag at 10 
/mi dCiardi et all 12005; Hais ch et all 12006). also does not sug- 
gest a large difference in extinction or bolometric luminosity be- 
tween the two stars. For the purposes of this paper, we assume 
that both objects are at a similar evolutionary state and embedded 
in an envelope. A third, faint comp onent in the SVS 20 system 
is separated from SVS 20 S by 073 dDuchene et aUl2007l) and is 
not detected in our spectra. 



Appendix B: Variable CO Ice Absorption from WL 6 

For WL 6, the depth of CO ice absorption decreased from 
r = 2.1 in the ISAAC spectrum to r = 0.51 in the CRIRES 
spectrum. From Gaussian fits to the absorption band, the central 
wavelength and FWHM remained similar, indicating no signifi- 
cant change in the ice composition. No other significant changes 
were detectable between the two observations. 

WL 6 is a poin t-s ource in K- and L-b an d AO 
imaging dRatzka et al.L 120051: iDuchene et all 120071) . 



lAlves de O liveira & Casali ( 2008) find variability of ~ 0.4 - 0.6 
mag in JHK monitoring of WL 6, which may be related to the 
variability in the CO ice absorption depth. If the disk is close 
to edge-on and occults the star only some of the time, then 
variability may be expected in the near-IR photometry and in 
the strength of the ice feature. 

Appendix C: Notes on CO Emission from Individual 
Targets 

In §3, we discussed the CO emission from individual stars in 
mostly generic terms. However, the interpretation of the narrow 
emission component from several objects is somewhat compli- 
cated. In the following subsections, we describe details of the 
narrow component from three objects, GSS 30, IRS 44 W, and 
CrA IRS 2. 

C.1. GSS 30 

Our CRIRES spectrum of GSS 30 covers 4.645-4.768 fjm and 
5.036-5.158 /im, with a large gap that excludes 12 CO lines with 
/' = 1 1 - 34 from our spectrum. Figure ICTI compares the 12 CO 
line profile for low-/ lines (/' < 10) to the scaled line profiles of 
12 CO with high-/ (35 < /' < 42), 13 CO, C ls O, and CO v = 2-1. 
The 12 CO low-/ lines are characterized by red- and blue-shifted 
emission with absorption at the cloud velocity. No absorption is 
detected in lines with / > 34. C I7 is detected in absorption but 
not in emission. 

The v = 2 — 1 lines can be described by the combination of 
a narrow, blueshifted Gaussian profile and a broader Gaussian 
profile centered at the systemic velocity. These two components 
roughly describe the emission profiles in lines of 13 CO, C ls O, 
12 CO with high-/, and I2 CO with low-/, though intervening 
absorption complicates the analysis. The broader component 10 is 
dominant in the high-/ I2 CO lines. The Gaussian profiles com- 
pared with the other lines in Figure IC.ll are kept in the same 
flux ratio as measured from the v = 2 — 1 lines. Excluding 
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Fig. C.l. A comparison of CO line profiles from GSS 30, as ex- 
tracted on-source. Two Gaussian profiles are fit to the emission 
in coadded v = 2 — 1 lines. The same two Gaussian profiles are 
then compared with the coadded lines of 12 CO v = 1 — with 
high-/, with low-/, 13 CO, and C 18 0. The two components in the 
v — 2 - 1 lines fit the other profiles reasonably well, although in 
the high-/ line the broader component dominates. Some narrow 
extended emission can also be seen as excess emission in low-/ 
I2 CO and I3 CO lines at -5 km s _1 . Both of these Gaussian com- 
ponents are considered narrow, optically-thick absorption within 
our crude classification scheme. 



the absorption, some minor differences can be seen between the 
Gaussian fits and the line emission. Very narrow emission in 
I2 CO and 13 CO transitions extends off-source in both the NE 
and SW within in the slit and is discussed in §3.3. 

To calculate the flux in each component, an equivalent width 
was measured over a specific velocity range on the blue and red 
side of each line. The equivalent width over this velocity range 
was then converted to a total flux based on Gaussian fits to the 
v -2-1 lines. This approach provides line fluxes measured with 
a consistent methodology for all lines of each isotope, despite 
that some model line profiles do not perfectly match the fit. 

Figure [8] shows the rotational diagram for narrow CO emis- 
sion from GSS 30, based on fits to the blue side of the emis- 
sion line. The flux ratio of 12 CO/ 13 CO lines is ~ 8, which in- 
dicates that the 12 CO lines are optically-thick. The flux ratio of 
I3 CO/C 18 lines of ~ 8 — 10 is simil ar to the abundance ratio of 
8.1 in the local ISM dWilsonL [l999h and suggests that both the 



10 Not to be confused with the broad component discussed in §3.1. 
Within our classification scheme, both components for GSS 30 are con- 
sidered narrow. 
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13 CO and C 18 emission are optically-thin. The measured tem- 
peratures for 13 CO and C ls O are 3 15 and 340 K. The 12 CO/ 13 CO 
line flux ratios yield a total CO column density log N(CO)= 18.6 
for b ~ 2.0 km s _1 and T = 340 K, assuming that the blueshifted 
12 CO emission is dominated by this component and not the 
broader emission component. At this column density, the ob- 
served 13 CO emission lines are nearly optically-thick, but not 
sufficiently enough to reduce the 13 CO/C 18 flux ratio. The 
fluxes in C 18 v = 1-0 and in 12 CO v = 2 — 1 lines do not suffer 
from any opacity effects until logA^(CO)> 19.5. At 340 K, the 
v = 2— 1 emission would not be detected, although a moderately 
higher temperature (~ 450 K) at this high column density would 
sufficiently populate the v = 2 level to detect v = 2 - 1 lines. The 
emission in high-/ lines of 12 CO and I3 CO is also not produced 
at such cool temperatures. The total emitting area is roughly (22 
AU) 2 for logJV(CO)= 18.6 + log(cosz). At a distance of 120 
pc, this spatial extent should be marginally resolvable. The ob- 
served CO emission is unresolved to a FWHM~ 0708 (10 AU at 
120 pc). The difference could be reconciled with a larger fe-value 
or a large viewing angle for our line-of-sight to the CO slab. On 
the red side of the line, the column density is difficult to assess 
because the non-detection of 13 CO emission is complicated by 
the red absorption. 

When analyzing spatially-extended CO e mission fr om GSS 
30 ob served with ISAAC (R = 10. 000). IPontoppidan etall 
(2002) found that the extraction region is much larger than the 
surface area implied by the line fluxes and column density. We 
confirm this problem with our CRIRE S data. To reconcile this 
discrepancy, Pontoppidan et al. (2002) suggested that the ex- 
tended emission is produced by material from the shock at the 
envelope/disk interface at 10-50 AU from the central star. This 
emission could then be reflected off the outflow cavity walls. 
This shock is expected to have a temperature plateau at ~ 500 
K due to reformation heating , with significant cooling in r ovi- 
brational CO and H 2 lines dNeufeld & Hollenbachi 1 1994 . In 
this scenario, the on-source and off-so urce line profiles should 
be similar, but IPontoppidan et al. (2002) was unable to test this 
prediction because the CO lines were unresolved with ISAAC. 
With the much higher resolution of CRIRES, we find that the 
off-source CO spectral line profiles are significantly narrower 
and have a higher peak-to-continuum ratio than the on-source 
CO line profiles. Although the extended nebulosity is bright in 
the K-band continuum, K-band pumping into v = 2 levels is 
ruled out by the lack of v = 2 - 1 emission detected off-source. 
An alternate but contrived explanation to reconcile the discrep- 
ancy in emission area versus extraction area is to invoke a very 
clumpy medium. Although the emission appears to be smoothly 
distributed within the slit, the spatial resolution may not be suf- 
ficient to detect many different clumps. 

C.2. IRS 44W 

The CO emission line profiles from IRS 44 W are similar to 
those from GSS 30. Figure [6] shows that I2 CO emission is seen 
on both sides of the CO absorption. However, the 13 CO and C ls O 
emission is seen only shortwards of the CO absorption. When 
the 13 CO emission line profile is scaled to the 12 CO profile (see 
inset in Fig. [6j>, a deficit in flux is seen in the 13 CO emission at 
< -30 km s _1 , at -5 km s (between the line peak and deep 
absorption), and on the entire red side of the line profile. 

The spatial profile of emission on the red wing of 12 CO lines 
also differs from the spatial profile of the blue wing of 12 CO 
lines and of the 13 CO lines. The emission on the red wing of 
12 CO lines is consistent with the location of continuum emis- 



sion from the secondary star. However, the 12 CO emission is lo- 
cated 0707 (8 AU at 120 pc) W of the continuum emission from 
the secondary star with a FWHM of ~ 0722 (26 AU at 120 pc). 
The blueshifted component is slightly stronger on the nights with 
poor seeing. The unresolved line equivalent widths are stronger 
in the ISAAC spectrum, which used a wider aperture and was 
obtained in worse seeing. These results both support the pres- 
ence of spatially-extended CO emission. The spectral and spatial 
information requires two physically distinct components for CO 
emission from IRS 44 W. 

For the blueshifted component, fits to the 13 CO and C 18 
line both yield temperatures of ~ 330 K. The 13 CO lines cannot 
be too optically thick, which places an upper limit on the column 
density of logA^(CO)< 19.4. The 12 CO lines are very optically- 
thick, with fluxes that should be considered upper limits, which 
places a lower limit on the column density of logA^(CO)> 18.9. 
At this column density, the approximate emitting area is ~ 6.5 
AU, less than that implied by the spatial extent of the emission. 

The red side of the line profile does not have reliable N(CO) 
and T because 13 CO and CO v = 2 — 1 were not detected. 

C.3. CrA IRS 2 

CrA I RS 2 is classified he re as an embedded object based on the 
SED dNutter et all 120051) . although insufficient evidence exists 
in the literature to confirm the presence of an envelope. 

The bright 12 CO and 13 CO emission lines from CrA IRS 2 
have non-Gaussian profiles with centroids that are redshifted by 
~ 2 km s _1 from the systemic velocity systemic and, when fit 
with a Gaussian profile, have FWHM of ~ 26 km s _1 . The v = 
1-0 lines also have a weak broad component that was discussed 
in §3.2. 

Line fluxes are calculated by fitting a median 13 CO line pro- 
file to every 12 CO and 13 CO line in the spectrum. The median 
13 CO line profile was calculated by coadding all 13 CO lines that 
are not affected by absorption. The resulting fits are typically 
good, except that high-/ I2 CO lines have broader peaks. The 
lack of any C I8 emission, with an upper limit of 15% of the 
flux in 13 CO lines, indicates that the 13 CO lines are not optically- 
thick. The 13 CO fluxes in the rotational diagram yield a best-fit 
temperature of T = 560 K. At this temperature, the 12 CO/ 13 CO 
line ratios indicate a column density log N(CO)~ 19.1, at about 
the limit where C ls O emission would be marginally detected. 
The total emitting area is (1.0 AU) 2 . 



Appendix D: Synthetic CO Line Fluxes 

For optically-thin gas, the temperature and number of CO 
molecules can be directly measured from an excitation dia- 
gram. Transitions become optically thick as the column density 
of the medium increases, thereby changing the line flux ratios 
and, as a consequence, direct temperature measurements. For 
CO emission from disks around CTTSs and Herbig AeBe stars, 
CO v = 1 - line flux ratios from disks deviate significantly 
from optically -thin branching ratios at a sing l e temperature (e.g . 
Britt ain et all 12003: Bla ke & BoogerA 12004 [Salvk et all [2007. 
2009). When only 12 CO lines are detected, such excitation dia- 
grams can be explained by emission produced in optically-thick 
gas or in gas with a large temperature gradient. The strengths of 
13 CO and CO v = 2 — 1 emission, relative to v = 1-0 emis- 
sion, suggests that the optically-thick interpretation is the best 
explanation for the curved rotational diagram obtained from our 
sample and from disks around CTTSs and Herbig AeBe stars. 
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We model the CO line fluxes by calculating the emission 
expected from an isothermal, ID slab of pure CO gas. The vi- 
brational excitation and rotational excitation of the CO gas are 
described b y the same temperatu re. The molecular data was ob- 
tained from Chandr a et al.ld 19961) . Throughout the layer, absorp- 
tion occurs in a Voigt profile with a Doppler broadening pa- 
rameter b, which includes thermal and turbulent broadening. 
Emission lines have a Gaussian profile with a FWHM of 1.82xfo. 
The density is assumed to be high enough so that collisions dom- 
inate the excitation, allowing the gas to maintain local thermal 
equilibrium. The layer is dust-free. The fraction of photons that 
escape from the total layer is calculated as a function of column 
density in each transition. 

The primary benefits of an analysis of CO line opacities are 
measuring accurate temperatures, approximate CO column den- 
sities, and rough surface areas for the emission regions. Within 
this paper, the Doppler parameter is set to b — 2.0 km s . A 
lower b parameter yields a faster increase in opacity, so resulting 
column densities would be smaller. The column densities in this 
simple ID slab model correspond to the amount of material at 
a specific temperature in our line of sight. For a slab model, the 
vertical column density would be log N = log N mo d + log cos i, 
where i is the incidence angle along our line of sight into the slab 
and Nmod is the column density of the model Characterizing this 
geometrical complication is beyond the scope of the simplistic 
approach adopted here. 

Figure ID. II shows the excitation diagram, normalized to 
(ij+\)A i °f ^ e ^(0) line, for a range of total CO column den- 
sities, T = 1000 K, and b — 1.0 km s . Each individual tran- 
sition starts to become optically thick at log Nj ~ 14 and is 
completely opaque at log Nj > 18. At 1000 K, the rotational 
population peaks at / = 13, so that transitions with J" ~ 13 
are the first to become optically thick. As a consequence, mid- 
/ line fluxes are weaker than expected, relative to low-/ and 
high-/ lines. The opacity in the R-branch transitions increases 
faster than the opacity in P-branch transitions because the lower 
levels of R-branch transitions are more populated than those 
in P-branch transitions. A more extreme example of the diver- 
gent P- and R-branch lines occ urs in the Orion BN/KL region 
ilGonzalez- Alfonso et al.Ll2002l) . 

For a slab with T = 1000 K, temperatures obtained from fits 
to synthetic optically-thick line ratios, either for the set of lines 
with /' < 30 or the combined set of lines with /' < 3 and 20 < 
/' < 30, are < 200 K lower than the input model temperatures. 
Temperature fits to only the high-/ (15 < / < 30) lines are 
tempting but can be misleading, with best fits as much as ~ 1000 
K higher than the input model temperature. Temperature fits to 
lines at / < 10 can yield temperatures of < 200 K and should be 
completely avoided. 
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